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INTRODUCTION 

Pulmonary  Function  Tests  (PFTs)  are  routinely  included  in  experimental  protocols  involving  hyper- 
and  hypobaric  excursions.  Unfortunately,  the  full  potential  of  those  tests  is  often  not  realized  due  to  the 
limited  capabilities  of  conventional  data  recording  technology.  A  potential  improvement  came  in  the  last 
decade  as  microcomputers  became  available  for  clinical  PFTs.  However,  microcomputer  use  in  the 
reseach  setting  has  been  limited  due  to  the  perceived  'black  box*  (i.e.  undocumented)  nature  of  many 
such  devices.  Nevertheless,  the  American  Thoracic  Society  (ATS)  stated  the  advantages  of  computer 
acquisition  of  pulmonary  function  data  (1)  as  follows: 

'(1)  Complete  automation  of  a  procedure  may  result  in  significantly  reduced  time  and  cost,  and  in 
increased  accuracy: 

(2)  Assurance  that  standardized  procedures  are  followed; 

(3)  Significant  reduction  in  major  measurement  errors; 

(4)  Storage  and  retrieval  of  information  quickly  and  efficiently; 

(5)  Implementation  of  automated  calibration  and  system  check  procedures  within  the  instrument; 

(6)  Standardized  and  consistent  interpretation  of  results.* 

A  recent  example  illustrates  the  power  In  these  statements.  Respiratory  flow  data  acquired  by  a 
computer  on  a  dive  to  450  m  revealed  oscillations  that  were  suspiciously  similar  to  the  frequency 
signature  of  HPNS  tremor  (2).  Had  a  conventional  mechanical  spirometer  been  used,  the  oscillations 
might  have  been  dismissed  as  mechanical  artifacts.  However,  the  accuracy  and  frequency  resolution  of 
the  computer  lent  credence  to  the  hypothesis  that  the  observed  oscillations  represented  a  previously 
unknown  action  of  the  High  Pressure  Nervous  Syndrome. 

The  ATS  suggests  comprehensive  quality  assurance  guidelines  for  the  use  of  digital  computers  in 
the  laboratory  (3).  Those  guidelines  were  refered  to  in  the  ATS  standards  for  spirometry  promulgated  in 
1987  (4).  According  to  the  standards,  users  of  PFT  software  should  possess  (and  hopefully  understand) 
source  code  listings  (computer  instructions  written  in  a  high  level  language  such  as  Basic,  C,  or 
Fortran).  Unfortunately,  commercial  software  is  typically  proprietary,  meaning  source  code  is  not 
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available.  Even  when  code  is  presented,  it  is  frequently  poorly  documented.  These  restrictions  are 
undesirable  in  the  clinical  setting,  but  are  intolerable  in  scientific  and  engineering  tasks. 

The  purpose  of  this  report  is  to  describe  the  hardware  installation  and  source  code  for  a  Pulmonary 
Function  data  acquisition  system  that  has  proven  useful  in  the  hyperbaric  research  laboratory.  It 
illustrates  how  computers  can  economically  collect,  analyze  and  display  data  from  a  variety  of 
pulmonary  function  tests  (Table  1).  [Readers  interested  in  computer  aided  measurement  of  dynamic, 
quasi-static  and  static  pulmonary  compliance  are  directed  to  Taytor  (5).] 
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METHODS 


Apparatus 

The  current  laboratory  uses  a  PC(AT)  done  (Model  Z-248,  Zenith  Data  Systems,  St.  Joseph,  Ml) 
although  any  IBM  compatible  system  would  suffice.  The  critical  component  is  the  analogue-to-digital 
(A/D)  convertor,  as  this  determines  both  the  maximal  sampling  rate  and  the  ease  of  software 
development.  The  MetraByte  DAS-16F  A/D  board  (Metrabyte  Corp.,  Taunton,  MA)  has  a  maximal 
sampling  rate  of  100  kHz,  can  take  single  ended  input  from  16  channels  (or  eight  differential  analogue 
inputs),  permits  both  foreground  and  background  data  collection,  and  is  fully  programmable  using  either 
commercially  available  or  laboratory-developed  software.  Row  was  measured  by  a  Hans  Rudolph 
pneumotachograph  (Model  4813,  Hans  Rudolph,  Inc.,  Kansas  City,  MO),  heated  to  37*Cand  coupled  to 
a  Validyne  differential  pressure  transducer  (MP45  ±  0.196  kPa  (±2  cmH^O),  Vaiidyne  Engineering, 
Northridge,  CA).  Mouth  pressure  was  obtained  using  a  second  Validyne  differential  pressure  transducer 
(MP45  ±  9.803  kPa  (±100  cmh^O))  connected  to  the  airway  between  the  mouth  and  the 
pneumotachograph.  Unfiltered  output  from  the  flow  and  pressure  transducers  was  amplified  using 
Validyne  CD18  and  CD19  carrier  demodulators,  respectively. 

The  frequency  content  of  flow  generated  during  the  Forced  Vital  Capacity  maneuver  extends  to  at 
least  12  Hz  (6).  Digital  sampling  theory  (7)  requires  sampling  to  occur  at  frequencies  at  least  double  that 
in  the  sampled  signal.  Accordingly,  data  from  both  inputs  was  sampled  at  50  Hz.  Sampling  at  even  100 
Hz  is  not  uncommon. 

[Validyne  is  now  manufacturing  an  interface  card  (UPC607)  which  serves  both  data  acquisition  and 
signal  conditioning  needs,  removing  the  need  for  separate  signal  demodulators]. 

Before  the  DAS-16F  board  can  be  used  It  must  first  be  correctly  installed  into  the  computer.  This 
installation  involves  both  hardware  and  software.  Physical  installation  requires  the  setting  of  four  D  I  P. 
switches,  and  then  the  correct  positioning  of  the  board.  The  user  must  correctly  position  the  switches 
for:  (a)  channel  configuration  (16  or  8  channel  use);  (b)  unipolar  (positive  only)  or  bipolar  (positive  and 
negative)  voltage  inputs;  (c)  direct  memory  access  level  (select  1  for  systems  with  hard-disk  and  3  for 
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(loppy  drive  systems);  and  (d)  gain,  which  operates  interactively  with  switch  (b)  to  set  the  input  voltage 
range.  Software  installation  is  performed  using  INSTALLEXE  (software  provided  with  the  board),  which 
sets  the  desired  base  address  for  the  board  within  your  system.  Once  installed  the  card  should  be 
calibrated  using  supplied  software  (CALF.EXE  or  CALBAS;  see  DAS-16F  Manual). 

Data  acquisition 

Appendix  1  contains  a  generic  subroutine  for  the  collection  of  data  in  either  foreground  mode  (all 
other  computer  functions  temporarily  suspended)  or  background  mode  (for  simultaneous  data  collection 
and  processing).  Data  acquisition  subroutines  are  presented  before  calibration  subroutines,  since  to 
calibrate  using  the  computer  you  must  first  be  able  to  collect  data.  All  computer  programs  were  written 
using  Microsoft  Fortran77,  and  are  presented  in  their  smallest  parts  (subroutines). 

The  information  below  provides  a  brief  overview  of  the  commands  issued  by  the  Fortran  software  to 
the  MetraByte  board.  Before  data  collection,  four  initializations  (operation  modes)  must  be  completed: 

(1)  Mode  0:  This  mode  is  used  to  initialize  the  A/D  board  by  first  testing  its 

physical  presence,  then  setting  an  input/output  (base)  address 
(variable  1),  setting  the  interrupt  level  (variable  2).  and  setting  the 
direct  memory  access  (D.M.A.)  level  (variable  3).  The  D  M  A.  level 
is  set  both  physically  on  the  board  and  within  the  software  (see 
DAS-16F  Manual). 

(2)  Mode  i:  This  mode  sets  the  lower  (variable  1)  and  upper  ranges  (variable  2) 

of  the  channels  to  be  scanned. 

(3)  Mode  17:  This  mode  sets  the  programmable  timer  rate  and  therefore 

determines  the  sampling  rate,  and  requires  two  variables:  (a)  divide 
data  for  counter  2  (variable  1),  and  (b)  divide  data  for  counter  1 
(variable  2).  These  counters  of  the  8254  counter  timer  produce  the 
output  pulse  rate  for  triggering  the  A/D  conversions.  The  sampling 
frequency  is  derived  from  the  following  equation: 
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Frequency  =  1,000,000/ [(variable  1 /number  of  channels)* 
variable  2]  However,  the  user  is  responsible  for  doing  this  division 
(a  point  not  in  the  manual),  so  set  the  size  of  variables  1  and  2 
with  this  in  mind.  For  example,  to  set  a  sampling  rate  of  50  Hz  for 
1  channel  variable  1  is  set  to  1000  and  variable  2  is  set  to  20.  To 
collect  2  channels  at  50  Hz,  leave  variable  2  at  20  and  set  variable 
1  to  500  (i.e.  1000/2). 

(4)  Mode  19:  This  last  initialization  mode  allows  for  the  use  of  an  external  trigger 

(positive  input  voltage)  to  initiate  the  start  of  A/D  conversions.  This 
is  an  extremely  useful  feature  which  permits  remote  operation  by 
either  the  experimenter  or  the  subject.  It  effectively  suspends  the 
computer  program  until  a  pre-determined  input  voltage  is  received. 

Three  variables  are  required  to  operate  this  mode. 

(a)  channel  number  for  the  analogue  input  voltage. 

(b)  the  magnitude  of  the  trigger  voltage  (bit  equivalent),  and 

(c)  the  slope  of  the  input  signal  (i.e.  positive  or  negative). 

Data  acquisition  is  performed  using  one  of  several  modes  (see  DAS-16F  Manual),  however  Modes  4 
(foreground)  and  5  (background)  are  sufficient  for  most  applications: 

(5)  Mode  4:  Performs  A/D  conversions  in  the  foreground  according  to  the 

channels  ranges  and  sampling  rates  determined  above.  Three 
variables  values  are  required:  (a)  the  number  of  conversions  (word 
count)  to  be  performed  (e  g.  for  2  channels  collecting  data  at  50 
Hz  for  10  seconds,  word  count  =  2  *  50  *  10);  (b)  an  array 
(previously  initialized  and  dimensioned  within  the  Fortran  code)  into 
which  the  input  shall  be  loaded,  and  the  position  within  the  array 
from  which  to  start  the  storage  (this  is  usually  zero  for  a  single 
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dimension  array  or  ‘[1,0]'  for  an  array  with  several  levels);  and  (c) 
the  source  of  the  trigger  to  start  the  A/D  conversions  (set  this  to  1 
to  use  the  previously  programmed  interval  timer,  Mode  17). 

(6)  Mode  5:  Performs  A/D  conversions  in  the  background  according  to  the 

channels  ranges  and  sampling  rates  determined  above.  Four  input 
variables  are  required:  (a)  the  number  of  conversions  (word  count) 
to  be  performed  (as  above);  (b)  the  segment  (16  bit)  of  the 
memory  to  receive  data  (you  need  to  seek  technical  help  here  as 
this  location  is  dependent  upon  the  memory  allocations  already 
assigned  within  your  system;  we  use  #5000  and  #6000  {Fortran}, 
or  &H5000  and  &H6000  {Basic});  (c)  the  source  of  the  trigger 
input  to  start  the  A/D  conversions  (as  above);  and  (d)  a  command 
telling  the  computer  to  carry  out  this  sequence  once  (single  cycle) 
or  to  operate  continuously  (until  stopped  using  Mode  7  -  see  DAS- 
16F  Manual). 

Calibration 

Calibration  of  the  pneumotachograph  must  be  performed  under  two  conditions;  zero  flow,  and 
using  a  standard  volume  (e.g.  3  liter)  passed  through  the  pneumotachograph  at  flow  rates  appropriate  to 
a  given  experiment.  Using  an  analysis  routine,  flow  is  integrated  with  respect  to  time,  providing  both 
volume  and  flow  calibration.  This  method  is  superior  to  electronic  integration,  since  the  latter  technique 
invariably  produces  u:i  electronic  drift  In  the  baseline,  making  data  interpretation  impossible. 

Mouth  pressure  is  often  of  interest,  especially  during  exercise  protocols.  Prior  to  calibration  the 
linearity  of  the  mouth  pressure  transducer  must  be  verified.  Subsequent  calibration  may  then  be 
performed  using  a  two  point  calibration  over  the  positive  range  of  anticipated  mouth  pressures.  From 
linear  regression,  the  pressure  data  yielded  both  the  zero  offset  and  slope  for  the  calibration.  These 
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values  were  saved  and  used  during  data  analysis.  The  calibration  subroutines  for  the  pneumotachograph 
and  mouth  pressure  transducer  are  contained  in  Appendix  2. 

Rest  and  exercise  respiratory  patterns 

Respiratory  pattern  data  (which  includes  VT,  VE,  and  respiratory  timing  variables)  were  derived  from 
the  same  routines  that  provided  lung  volume  data.  The  difference  lay  in  the  breathing  pattern  used  by 
the  subject,  rather  than  within  the  computer  program.  The  example  provided  In  Appendix  3  is  from  a 
program  written  to  analyze  data  collected  at  450  m  (changes  in  ambient  pressure  are  automatically 
taken  into  account  by  these  procedures).  Subjects  exercised  for  25  min  while  exposed  to  a  moderate 
external  airway  resistance.  Appendix  8  contains  the  entire  program.  This  program  also  included  the 
measurement  of  mouth  pressure  (including  the  determination  of  peak-to-peak  values  and  inspiratory  to 
expiratory  ratios)  and  the  integration  of  mouth  pressure  into  trajectory  data  for  a  flow-volume-pressure 
three-dimensional  surface  (8). 

Lung  volume  measures 

The  measurement  of  lung  volumes  (exc  :ding  RV)  relies  more  upon  the  skills  of  the  experimenter 
than  the  art  of  the  programm.  r  ''ubjects  must  be  trained  to  provide  reproducible  VC  maneuvers  through 
the  attainment  of  consistent  RV  and  TLC  end  points.  Once  this  is  accomplished  the  experimenter  can, 
with  a  single  maneuver,  obtain  VC,  ERV,  1C,  and  IRV  (Figure  1).  Following  an  exhalation  to  RV,  the 
subject  slowly  inspired  to  maximum  lung  capacity  (TLC),  followed  by  an  expiration  to  RV  and  a  normal 
slow  inspiration.  The  subject  then  took  normal  breaths  (giving  VT,  ERV,  and  by  subtraction,  1C).  From 
this  point  the  subject  inspired  to  TLC  again.  By  subtracting  from  1C  the  VT  obtained  using  Appendix  3 
(since  the  VT  in  this  maneuver  is  often  falsely  large),  IRV  is  derived. 

Subroutines  necessary  for  measuring  the  above  lung  volumes  are  contained  in  Appendix  4.  Note 
that  the  programming  code  is  the  same  for  each  lung  volume;  the  changes  being  only  in  the  maneuver 
performed  by  the  subject.  These  subroutines,  and  those  included  below  (except  the  flow-volume 
subroutines),  perform  analysis  on  complete  respiratory  cycles  only.  This  is  essential,  since  during 
exercise  subjects  do  not  commence  inspiration  at  the  same  time  as  the  computer  starts  to  collect  data, 
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and  rarely  terminate  exhalation  immediately  prior  to  the  end  of  data  collection.  Therefore,  within  the  initial 
analysis  subroutine  (Appendix  3)  is  a  'paragraph'  which  searches  the  flow  data  to  find  the  first 
commencement  of  inspiration.  This  is  found  when  the  flow  data  crosses  the  predetermined  basal  flow 
threshold  in  the  inspiratory  direction.  The  subroutine  then  counts  the  number  of  inspirations  (counter  i) 
and  expirations  (counter  E).  To  ensure  these  are  always  equal,  the  last  inspiration  is  omitted,  and  the 
output  is  always  written  for  1-1  inspirations  and  expirations. 
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FIGURE  1 


PRIMARY  SUBDIVISIONS  OF 
LUNG  VOLUME 


SPECIAL  DIVISIONS  FOR 
PULMONARY  FUNCTION  TESTS 


Standard  lung  volumes  and  capacities. 


Maximal  voluntary  ventHation 


The  MW  maneuver  (Appendix  5)  is  also  a  hybrid  of  the  breathing  pattern  subroutines  (Appendix  3). 
For  the  MW,  all  the  VTs  are  summed  then  divided  by  the  time  taken  to  complete  the  maneuver.  This 
time  is  controlled  by  the  experimenter  so  that  It  equals  the  duration  of  interest  (i.e.  12  or  15  seconds).  As 
with  the  above  subroutines,  only  complete  cycles  are  analyzed,  so  the  time  to  complete  the  maneuver 
(as  determined  by  the  software)  w8l  often  not  equal  the  time  used  during  the  experiment.  This  difference 
is  minimal,  and  results  in  greater  precision  in  the  determination  of  MW  than  is  possible  using  a 
spirometer. 

Row-volume  loops  and  derived  variables 

For  this  maneuver  (Appendix  6)  subject  training  was  again  vital  to  ensure  that  subjects  obtained 
true  maximal  flows  and  reproducible  end  points  at  RV  and  TLC.  The  amplifier  gain  setting  was  adjusted 
so  that  voltages  during  peak  flows  were  as  dose  as  possible  to  the  _+5  Volt  input  capacity  of  the  A/D 
board. 

Since  data  sampling  was  at  50  Hz,  flow  and  instantaneous  lung  volumes  were  derived  at  20  ms 
intervals  over  the  entire  volume  excursion.  These  data  were  saved  as  an  ASCII  file  for  importation  into 
suitable  graphics  packages  (e.g.  Lotus123  (Figure  2)  or  Sigmaplot  [see:  Taylor  1990]). 

The  flow  data  obtained  from  these  maneuvers  was  also  used  to  obtain  the  FEV,  0  for  each 
expiration.  The  ASCII  flow  data  was  searched  manually  to  detect  PEFR  and  MEFR  for  various 
percentages  of  the  VC.  While  algorithms  for  these  latter  computations  may  be  written  and  incorporated 
into  the  analysis  subroutines,  it  is  often  easier  to  take  these  values  by  hand  directly  from  the  ASCII  files. 

If  noise  occured  in  the  flow  signal  at  the  volume  of  interest,  computer  analysis  could  be  confounded. 
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APPENDIX  1 


Data  acquisition  subroutines. 


ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 

c  GENERIC  DATA  COLLECTION  PROGRAM 

c 

c  Program  uses  Metrabyte  A/D  conversion  board,  collects  data 
c  at  50  Hz  from  channels  3  and  4.  Collection  occurs  in  the  background 
c  (Mode  5)  and  takes  place  in  1  minute  segments.  Between  each  segment 
c  the  program  looks  at  channel  3  for  analog  input  from  an  external 
c  trigger.  This  input  terminates  the  data  collection.  Data  is  continuously 
c  written  to  disk  in  ASCII  code  between  successive  samplings, 
c 

c  NIGEL  A  S.  TAYLOR  4th  MAY  1990 

c  DIVER  UFE  SUPPORT  EQUIPMENT  PROGRAM,  NAVAL  MEDICAL  RESEARCH  INSTITUTE 
c 

c  LATEST  UPDATE:  5/24/90 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 

$STORAGE:2 

INTEGER  PARAM,IM0DE,0FFADR,DASH16,RC0DE,J,I,L,M,N,0,ZZ 
INTEGER  DATA1  ,DATA2,LMODE,LCODE,LP ARAM, JCODE.JMODE, COUNT 
INTEGER  ANSI  .ANS2,TIME,CHECK,KCODE,KMODE,VAR,ELEM,TCOUNT 
DIMENSION  DATA1  (3,3001  ),PARAM(4),DATA2(3, 3001) 

DIMENSION  LPARAM(5),ELEM(5),VAR(5) 

CHARACTER*20  FILENAME 
CHARACTER  BELL 
BELL  =  CHAR  (7) 
c 

c  The  bell  is  used  to  produce  a  'beep'  from  the  PC.  This  prompts 
c  the  experimenter  for  input.  See  the  WRITE  statements  below, 
c 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 

OPENH 

100  WRITE(*,l) 

I  FORMAT  (/////////// 

•  *  ENTER  NAME  OF  OUTPUT  RLE  FOR  DATA :'/////////) 

READ(MI)  FILENAME 

II  FORMAT(A20) 

OPEN(3,FILE  =  FILENAME, STATUS = 'NEW,  FORM = 'FORMATTED') 

WRITE(*,7000) 

7000  FORMAT^/  ENTER  EXPERIMENTAL  DURATION  IN  MINUTES  (13):'//) 

READ(*,7001)  TIME 

7001  FORMAT(l6) 

TCOUNT =0 
CHECK =0 

zz*o 

c 
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c  Initialization  of  the  MetraByte  A/D  board:  see  text  for  description 
c 

IMODE  =  0 
PARAM(1 )  «  #300 
PARAM(2)  *  2 
PARAM(3)  *  3 

RCODE  =  DASH16(IMODE,  PARAM) 
c 

c  Setting  the  channel  scannii.g  limits  for  3  channels,  starting  at  2  and 
c  ending  at  4. 
c 

IMODE  =  t 
PARAMO)  =  2 
PARAM(2)  -  4 

RCODE  =  DASH16(IMODE,  PARAM) 

c 

c  Setting  the  sampling  rate  for  three  channels  at  50  Hz 
c  (i.e.  1.000,000/[{333=  1000/3}*20)) 


IMODE  =  17 
PARAMO)  =  333 
PARAM(2)  =  20 

RCODE  =  DASH16(IMODE,  PARAM) 

WRITE(M2) 

12  FORMAT (///■  PRESS  THE  TRIGGER  TO1/ 

*  •  ******  ACTIVATE  DATA  COLLECTION  ******•/////) 


c 


c  Using  an  analog  trigger  on  channel  2  to  activate  data  collection, 
c  Program  pauses  until  an  positive  voltage  greater  than  300  (binary 
c  units)  is  detected.  This  is  a  very  useful  Mode  that  can  be  used  at 
c  any  point  in  the  program.  Unlike  large  machines,  the  PC  is  unable 
c  to  check  hardware  status  while  completing  other  tasks,  so  this  Mode 
c  can  also  be  used  to  terminate  operations  as  well  as  to  initiate. 


c 


c 
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IMODE  =  19 
PARAM(l)  =  2 
PARAM(2)  =  300 
PARAM(3)  =  0 

RCODE  =  DASH16/IMODE,  PARAM) 

WRITE(M41)  BELL 
WRITE/*, 13) 

FORMAT(//// 

. ***** . 7 

•••COLLECTION  STARTED  •••’/ 

. **** . ***///) 


c 

c  Background  data  collection  using  Mode  5.  See  test  for  description, 
c  Note  that  the  Variable/Arrays  for  this  CALL  have  been  altered  from 
c  previous  modes.  This  is  because  they  need  to  be  unique  to  this  CALL 
c  since  it  will  be  operating  in  the  background  and  will  be  confused  if 
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c  these  values  are  altered  during  operation, 
c 

LMODE  =  5 
LPARAM(l)  =  9000 
LPARAM(2)  =  #5000 
LPARAM(3)  =  1 
LPARAM(4)  =  0 

LCODE  =  DASH16(LMODE,  LPARAM) 

ZZ  =  ZZ  +  1 
c 

7777  CONTINUE 
c 

c 

c  Mode  8  is  used  to  monitor  the  status  of  the  background  data  collection, 
c  In  this  example  we  wish  1  full  minute  of  data  collection  to  be  transferred 
c  into  an  array  from  memory  (3  chans  *  50  Hz  *  60  sec  =  9000).  Mode  8  keeps 
c  the  program  sitting  till  the  Word  count'  equals  9000.  We  could  use  this 
c  Mode  to  transfer  data  in  smaller  segments  for  real-time  processing.  This  is 
c  the  real  advantage  that  Mode  5  provides, 
c 

3330  JMODE  =  8 

JCODE  =  DASH  16( JMODE,  VAR) 

IF  (VAR(3).GE.9000)  GOTO  3001 
GOTO  3330 

3001  WRITEC3900)  ZZ 

3900  FORMAT('+  MINUTE  ,12,'  COMPLETE:  ) 
c 

c  Here  Mode  4  is  used  to  monitor  channel  2  for  a  positive  input  from 
c  the  trigger  (i.e.  greater  than  binary  value  of  500).  Mode  1  is  not 
c  re-defined,  instead  all  3  channels  are  monitored.  However,  the  Word 
c  count'  is  set  to  1  to  remove  the  input  from  the  other  channels, 
c  Note:  this  only  works  because  the  trigger  is  connected  to  channel  2. 
c  If  trigger  input  is  detected  the  program  jumps  to  Statement  4000  to 
c  save  the  data,  then  to  Statement  8  to  terminate  data  collection, 
c 

IMODE  =  4 
PARAM(l)  =  1 

PARAM(2)  *=  OFFADR(CHECK) 

PARAM(3)  =  1 

RCODE  -  DASH  16(IMODE,  PARAM) 

IF  (CHECK.GT-500)  GOTO  4000 
IF  (TCOUNT.GE.TIME)  GOTO  4000 


c  Mode  5  is  activated  again  immediately  after  the  trigger  check.  This 
c  results  in  minimal  data  loss  (milliseconds), 
c 

LMODE -  5 
LPARAM(l)  -  9000 
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LPARAM(2)  =  #6000 
LPARAM(3)  =  1 
LPARAM(4)  =  0 

LCODE  =  DASH16(LMODE,  LP ARAM) 

ZZ  =  ZZ  +  1 
c 

c  Mode  9  is  then  used  to  transfer  data  from  memory  to  a  previously 
c  initialized  and  dimensioned  array.  This  occurs  in  the  foreground 
c  while  the  data  acquisition  subroutines  are  operating  in  background, 
c  Five  input  Variables  are  required  for  Mode  9: 
c  (1)  *word  count',  9000  =  3  chans  *  50  Hz  *  60  seconds 
c  (2)  memory  address  as  per  Mode  5  declaration 
c  (3)  the  conversion  number  to  start  the  transfer  from 
c  (4)  the  array  name  and  location  within  the  array  (e.g.  point  0  of  level  1) 
c  (5)  channel  array  data:  use  0  if  not  required, 
c 

4000  KMODE  =  9 

ELEM(l)  =  9000 
ELEM(2)  =  #5000 
ELEM(3)  =  0 

ELEM(4)  =  OFFADR(DATA1(1,0)) 

ELEM(5)  =  0 

KCODE  =  DASH16(KMODE,  ELEM) 
c 

c  Now  while  Mode  5  is  still  operating  the  data  are  written  to  a 
c  file  on  the  disk  in  ASCII  code.  This  file  is  not  closed,  but  is 
c  continually  added  to  with  data  from  each  usuccessive  minute, 
c 

DO  3002  I  =  1,(9000/3) 

WRITE(3,4949)  (DATAl(J,iy  =  23) 

3002  CONTINUE 
4949  FORMAT(2I6) 

c 

TCOUNT  =  TCOUNT  +  1 
IF  (CHECK.GT 300)  GOTO  8 
IF  (TCOUNT.GE.TTME)  GOTO  8 
c 

3331  JMODE  =  8 

JCODE  =  DASH16(JMODE,  VAR) 

IF  ( VAR(3)  .GE.9000)  GOTO  3004 
GOTO  3331 ' 

3004  WRITEC3900)  ZZ 
c 

IMODE  =  4 
PARAM(l)  =  1 

PARAM(2)  =  OFFADR(CHECK) 

PARAM(3)  =  1 

RCODE  =  DASH  16(IMODE,  PARAM) 
c 

IF  (CHECK.GT 300)  GOTO  4001 
IF  (TCOUNT.GE.TIME)  GOTO  4001 


c 
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c 

4001 


c 


3005 

c 


c 


8 

141 

14 


20 

c 


c 


LMODE  =  5 
LPARAM(l)  =  9000 
LPARAM(2)  =  #5000 
LPARAM(3)  =  1 
LPARAM(4)  =  0 

LCODE  =  DASH16(LMODE,  LPARAM) 
ZZ  =  ZZ  +  1 


KMODE  =  9 
ELEM(l)  =  9000 
ELEM(2)  =  # 6000 
ELEM(3)  =  0 

ELEM(4)  =  OFFADR(DATA2(1,0)) 
ELEM(5)  =  0 

KCODE  =  DASH  16(KMODE,  ELEM) 

DO  3005  1  =  1,(9000/3) 

WRITE(3,4949)  (DATA2(J,I)rJ  =  23) 
CONTINUE 

TCOUNT  =  TCOUNT  +  1 
DF  (CHECK.GT 300)  GOTO  8 
IF  (TCOUNT.GE.TTME)  GOTO  8 

GOTO  7777 


WRITE(*,141)  BELL,BELL,BELL,BELL,BELL 
FORMAT('  \A1,'  ’Ah'  Ah'  'Ah'  \A1) 
WRITER, 14) 

FORMAT(// 


*•*  COLLECTION  FINISHED 


FORMAT(Il) 


'*•••  j 

■*•*'/) 


CLOSE(3) 


WRITE(*3101)  FILENAME 
2101  FORM.\T(/ ///////// 

*  '  DATA  FROM:  ',A2Q,'  SAVED.'/// 

*  . 7 

*  '  RECORD  EXPERIMENT  TIME!!'/ 

*  . 7/ 

*  '  WHAT  DO  YOU  WISH  TO  DO  NEXT?'/ 

*  '  1  RETURN  =  COLLECT  MORE  DATA'/ 

*  '  9  RETURN  =  TERMINATE  PROGRAM’/) 

READ(\20)  ANS2 

IF  (ANS2.EQ.1)  GOTO  100 


17 


IF  (ANS2.EQ.9)  CONTINUE 
c 

9999  WRITER, 9191) 

9191  FORMAT  (///////////// 


STOP 

END 


GOOD  BYE 


*7 

7 

•'/////) 


APPENDIX  2 


Calibration  subroutines. 


PART  I:  Collecting  the  calibration  data. 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccaxccccccccmxxcccccccccccccccccccc 

c 

c  EXAMPLE  CALIBRATION  Program: 

c 

c  Refer  to  the  GENERIC  DATA  COLLECTION  program  for  explanations 
c  of  various  programming  CALL  routines  and  functions, 
c 

c  CALIBRATION  FOR  MOUTH  PRESSURE,  END-TIDAL  C02  AND 
c  PNEUMOTACHOGRAPH 

c 

c  SET-UP  FOR  THE  DASH  16  A/D  BOARD  (MetraByte) 
c  Data  collection  started  from  input  to  channel  two. 
c  Sampling  at  50Hz. 

c  Four  channels  monitored:  2  =  trigger,  3= mouth  pressure 
c  4 = flow,  5  =  end-tidal  C02. 

c 

c  NIGEL  A.S.  TAYLOR  May  1990 

c  Diving  Medicine,  Naval  Medical  Research  Institute, 
c 

c  LATEST  UPDATE:  21  MAY  1990 

cxxxxcccccccccccccccccccccccccccccccccaxccccccccccccccccccccccccccccccccccccccc 

$STORAGE:2 

INTEGER  ZERO,GAIN,PEAKMODE,ZCHAN,HIGH,MM,M,B,D,E,F,G,HrI,K,IMODE 
INTEGER  ZER,HIG, COUNT, GAIN2,BASE,BAS,ICHAN, BIG, BOG.PARAM, SAVE 
INTEGER  MEANNOFLOW.LITER, NOFLOW, GAIN3,ANS,DASH16,OFFADR,RCODE 
REAL  MEANZERO.MEANHIGH, PRESS  UR  E,REGCOEFF, CONST  ANT 
REAL  V,W,X,Y,Z 

REAL  RZERO,RHIGHAC,RBASE,RBIG,PRESS2 
CHARACTER  *40  FILENAME, FILE2,FILE3,F1LE4 
CHARACTER  BELL 

DIMENSION  ICHAN(2),ZERO(ll),HIGH(ll),RZERO(ll),RHIGH(ll) 

DIMENSION  BASE(11),RBASE(11),BIG(11),RBIG(11),PARAM(4) 

DIMENSION  NOFLOW(101)rUTER(2000) 

cccccccccccccccccccccccccceccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

OPEN(') 

BELL  =  CHAR(7) 

100  WRITE(*,1) 

1  FORMAT  (/////////////// 

•  '  PRESSURE,  PetC02  AND  VOLUME  CALIBRATION  Program.'////// 

•  '  NOTE:  YOUR  FIRST  CALIBRATION  MUST  BE  THE'/ 

•  '  TRANSDUCER  MEASURING  MOUTH  PRESSURE.'// 

•  '  ENTER  1  [RETURN]  TO  COMMENCE  CALIBRATION///) 

READC.20)  ANS 


20 

c 


IF  (ANS.EQ.1)  CONTINUE 
FORMAT(Il) 


cccca:ccccccccccccccccca:ccccccaxx:ccccccccxxccccccaxccccctxccccaxccxxcaxccccccccc£caxcccccccc 
c  Pressure  calibration 
c 

10  IMODE  =  0 

PARAM(l)  =  #300 
PARAM(2)  =  2 
PARAM(3)  =  3 

RCODE  =  DASH  16( IMODE,  PARAM) 
c 

IMODE  =  1 
PARAM(l)  =  3 
PARAM(2)  =  3 

RCODE  =  DASH16(IMODE,  PARAM) 
c 

IMODE  =  17 
PARAM(l)  =  1000 
PARAM(2)  =  20 

RCODE  =  DASH16(IMODE,  PARAM) 
c 

WRITEC3) 

3  FORMAT(///'  SET  THE  PRESSURE  TRANSDUCER  AT  ZERO  PRESSURE.’/ 

*  IS  TRIGGER  CONNECTED  TO  CHANNEL  TWO?'/ 

•  '  IS  TRANSDUCER  CONNECTED  TO  CHANNEL  THREE?’///// 

.  .  .....  pRESS  TRIGGER  TO  COMMENCE  *•*•*’//) 

c 

IMODE  =  19 
PARAM(l)  =  2 
PARAM(2)  =  300 
PARAM(3)  =  0 

RCODE  =  DASH  16(IMODE,  PARAM) 
c 

WRITE(*31) 

31  FORMAT(/’ . COLLECTION  HAS  BEGUN . ’/) 

WRITE(*,4000)  BELL 
c 

IMODE  =  4 
PARAM(l)  =  12 

PARAM(2)  =  OFFADR(ZERO(0)) 

PARAM(3)  =  1 

RCODE  =  DASH16(IMODE,  PARAM) 

c 

c  Summing  10  zero  pressure  points  to  obtain  the  mean  (below), 
c 

A  =  0.0 
DO  311  J  =  2,11 
ZER  =  ZERO(J) 

RZERO(J)  =  FLOAT(ZER) 

A  =  A  +  RZERO(J) 
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311  CONTINUE 
c 

WR1TE(*,4) 

4  FORMAT(/' . CALIBRATION  VALUES  COLLECTED . '/) 

WRTTE(*,4000)  BELL 

c 

c  Calculating  mean  zero  pressure, 
c 

MEANZERO  =  A  /  10.0 

WRITE(*,41)  (RZERO(J),J  =  2,11), MEANZERO 

41  FORMAT(2(5(2X4:9.2)///V  MEAN  ZERO  PRESSURE  =  ',F9 2/f) 

WRITE(*,42) 

42  FORMAT('  CHECK  BASAL  VALUES  ARE  ALL  ABOUT  EQUAL.  ) 
c 

WRITE(*,43) 

43  FORMAT(//'  WHAT  DO  YOU  WISH  TO  DO  NEXT?'/ 

*  '  1  RETURN  =  CONTINUE  TO  NEXT  STAGE,'/ 

*  '  6  RETURN  =  REPEAT  &  ERASE  BASAL  CALIBRATION,'/ 

*  '  9  RETURN  =  END  PROGRAM.'/) 

READ(*,20)  B 

IF  (B.EQ.l)  CONTINUE 
IF  (B.EQ.6)  GOTO  10 
IF  (B.EQ.9)  GOTO  9 
c 

50  WRITE(*,5) 

5  FORMAT(///'  APPLY  SUITABLY  LARGE  PRESSURE  TO  THE  TRANSDUCER.'/ 

*  THIS  PRESSURE  SHOULD  BE  GREATER  THAN  THE  PRESSURE'/ 

*  '  YOU  ANTICIPATE  DURING  THE  EXPERIMENT.'/////) 

c 

IMODE  =  0 
PARAM(l)  =  #300 
PARAM(2)  =  2 
PARAM(3)  =  3 

RCODE  =  DASH  16(IMODE,  PARAM) 
c 

IMODE  =  1 
PARAM(l)  =  3 
PARAM(2)  =  3 

RCODE  =  DASH16(IMODE,  PARAM) 
c 

MODE  *  17 
PARAM(l)  =  1000 
PARAM(2)  =  20 

RCODE  =  DASH16(IMODE,  PARAM) 
c 

WRITE(*,51) 

51  FORMATO  *•**  PRESS  TRIGGER  TO  COLLECT  CALIBRATION  •*••'//) 
c 

IMODE  =  19 
PARAM(l)  =  2 
PARAM(2)  =  300 
PARAM(3)  =  0 
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RCODE  =  DASH  16(IMODE,  PARAM) 

WRITER) 

WRITE(*,4000)  BELL 
c 

IMODE  =  4 
PARAM(l)  =  12 

PARAM  (2)  =  OFFADR(HIGH(0)) 

PARAM(3)  =  1 

RCODE  =  DASH16(IMODE,  PARAM) 
c 

C  =  0.0 
WRITE(*,4) 

WRITE(  *  ,4000)  BELL 
c 

c  Summing  10  peak  pressures, 
c 

DO  511  K=2,ll 
HIG  =  HIGH(K) 

RHIGH(K)  =  FLOAT(HIG) 

C  =  C  +  RHIGH(K) 

511  CONTINUE 
c 

c  Calculating  mean  peak  pressure. 

c 

MEANHIGH  =  C  /  10.0 

WRITE(*,52)  (RHIGH(K),K = 2,11),MEANHIGH 

52  FORMAT(2(5(2X,F9.2)///),'  MEAN  HIGH  PRESSURE  =  VF9.2//) 

WRITE(*,53)  MEANZERO 

53  FORMATC  MEAN  ZERO  PRESSURE  =  \F9.2//) 
c 

WRITE(*,54) 

54  FORMAT(//'  WHAT  DO  YOU  WISH  TO  DO  NEXT?'/ 

*  '  1  RETURN  =  LINEAR  REGRESSION  ANALYSIS,'/ 

*  ’  6  RETURN  =  REPEAT  POOR  HIGH  PRESSURE  CALIBRATION,'/ 

*  ’  9  RETURN  =  END  Program.'//) 

READ(*,20)  D 

IF  (D.EQ.l)  CONTINUE 
IF  (D.EQ.6)  GOTO  50 
IF  (D.EQ.9)  GOTO  9 
c 

WRITER, 6) 

6  FORMAT(//'  PLEASE  TYPE  THE  PRESSURE  OF  THE  SECOND  CALIBRATION'/ 

*  '  USING  cm.H20  (e.g.  xxxxjcx  [F7.2]) '//) 

READ(*,61)  PRESSURE 

61  FORMAT(F7.2) 


c  Linear  regression  analysis  to  obtain  coefficients  for 
c  equation:  y  =  A  +  Bx.  Where:  A  =  y  intercept, 
c  B  =  slope, 

c 
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W  =  2  *  PRESSURE  *  MEANHIGH 
X  =  PRESSURE  *  (MEAN ZERO  +  MEANHIGH) 

Y  =  2  *  PRESSURE  *  PRESSURE 
Z  =  PRESSURE  *  PRESSURE 
REGCOEFF  =  (W  -  X)  /  (Y  -  Z) 

V  =  MEANZERO  +  MEANHIGH 
CONSTANT  =  (V  -  REGCOEFF  *  PRESSURE)  /  2 

c 

WRITE(*,7)  REGCOEFF, CONSTANT 

7  FORMAT(////’  LINEAR  REGRESSION  ANALYSIS  HAS  BEEN  PERFORMED.'/ 

*  '  REGRESSION  COEFFICIENT  B  (cmH20)  =  \F 9.2 / 

*  '  CONSTANT  A  (cmH20)  =  ,F9.2//) 

c 

71  WRTTE(*,72) 

72  FORMAT(///'  TO  CALCULATE  ACTUAL  PRESSURES  USE  THE  EQUATION:'/ 

*  x  =  (computer  raw  unit  -  A)  /  B.1/// 

*  1  ENTER  NAME  OF  PRESSURE  CAL  FILE:  //) 

WRITER  *,4000)  BELL 

READ(*,73)  FILENAME 

73  FORMAT(A40) 
c 

OPEN(3,FILE  =  FILENAMESTATUS  =  NEW) 

WRITE(3,8)  REGCOEFF, CONSTANT 

8  FORMAT(2X,F9.2,2X,F9.2) 

CLOSE(3) 

c 

821  WRITE(*,83) 

83  FORMAT(///'  WHAT  DO  YOU  WISH  TO  DO  NOW?'/ 

*  '  1  RETURN  *  CALIBRATE  THE  MASS  SPEC,'/ 

*  6  RETURN  =  CALIBRATE  PNEUMOTACHOGRAPH,'/ 

*  '  9  RETURN  =  END  Program.'//) 

READ(*,20)E 

IF  (E.EQ.l)  CONTINUE 
IF  (E.EQ.6)  GOTO  86 
IF  (E.EQ.9)  GOTO  9 

ccccccccccccccccccccccccccccmxcccccccccccccccccxxccccxcccccccccccccccccccccccc 
c  End-tidal  C02  calibration 

c 

c  This  calibration  is  exactly  the  same  as  that  used  for  the  mouth 
c  pressure  calibration  subroutines  (computers  are  not  smart  enough 
c  to  know  the  differOice!!).  Therefore  no  comments  are  provided, 
c 

84  IMODE  =  0 
PARAM(l)  =  #300 
PARAM(2)  =  2 
PARAM(3)  =  3 

RCODE  =  DASH16(IMODE,  PARAM) 
c 

IMODE  =  1 
PARAM(l)  =  5 
PARAM(2)  =  5 

RCODE  -  DASH  16(IMODE,  PARAM) 
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IMODE  =  17 
PARAM(l)  =  1000 
PARAM(2)  =  20 

RCODE  =  DASH  16(IMODE,  PARAM) 
c 

WRITER, 8400) 

8400  FORMAT(///'  SET  THE  MASS  SPEC  OUTPUT  TO  ZERO.'/ 

*  1  CHECK  TRANSDUCER  CONNECTED  TO  CHANNEL  5.'/ 

.  .  PRESS  TRIGGER  TO  COMMENCE  •••**■//) 

c 

IMODE =19 
PARAM(1)  =  2 
PARAM(2)  =  300 
PARAM(3)=0 

RCODE  =  DASH16(IMODE,  PARAM) 
c 

WRITEC31) 

WRITER, 4000)  BELL 

c 

IMODE =4 
PARAM(1)  =  12 

PARAM(2) = OFFADR(BASE(0)) 

PARAM(3)  =  1 

RCODE  =  DASH  16(IMODE,  PARAM) 

c 

A  =  0.00 
DO  841  J  =  2,ll 
BAS  =  BASE(J) 

RBASE(J)  =  FLOAT(BAS) 

A  =  A  +  RBASE(J) 

841  CONTINUE 
c 

WRITE(*,4) 

WRITE(  *  ,4000)  BELL 
MEANZERO  =  A  /  10.00 
WRITE(*,41)  (RBASE(jy  =  2,11), MEANZERO 
WRITE(*,42) 
c 

WRITE(*,43) 

READ(*,20)  G 
IF  (G.EQ.l)  CONTINUE 
IF  (G.EQ.6)  GOTO  84 
IF  (G.EQ.9)  GOTO  9 
c 

WRITE(*^111) 

5111  FORMAT(//’  APPLY  SUITABLY  HIGH  C02  OUTPUT  FROM  MASS  SPEC.  //) 
c 

8420  IMODE  =  0 

PARAM(l)  =  #300 
PARAM(2)  =  2 
PARAM(3)  =  3 
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RCODE  =  DASH  16(IMODE,  PARAM) 
c 

IMODE  -  1 
PARAM(l)  -  5 
PARAM(2)  =  5 

RCODE  =  DASH16{IMODE,  PARAM) 
c 

IMODE  =  17 
PARAM(l)  =  1000 
PARAM(2)  =  20 

RCODE  =  DASH16(IMODE,  PARAM) 
c 

WRITEC^l) 

c 

IMODE  =19 
PARAM(1)  =  2 
PARAM(2)  =  300 
PARAM(3)=0 

RCODE  =  DASH16(IMODE,  PARAM) 
c 

WRITE(*,31) 

WRITE(*,4000)  BELL 
c 

IMODE=4 
PARAM(1)  =  12 

PARAM  (2)  =  OFFADR(BIG(0)) 

PARAM(3)  =  1 

RCODE  =  DASH16(IMODE,  PARAM) 
c 

WRITEC,4) 

WRITE( *,4000)  BELL 

c 

C  =  0.0 

DO  843  K=2,ll 
BOG  =  BIG(K) 

RBIG(K)  =  FLOAT(BOG) 

C  =  C  +  RBIG(K) 

843  CONTINUE 
c 

MEANHIGH  =  C  /  10.0 
WRITE(*,52)  (RBIG(K),K  =  2,11), MEANHIGH 
WRITE(*43)  MEAN  ZERO 
c 

WRITE(V54) 

READC.20)  H 
IF  (H.EQ.l)  CONTINUE 
IF  (H.EQ.6)  GOTO  8420 
IF  (H.EQ.9)  GOTO  9 
c 

WRITE(*,6111) 

6111  FORMAT(///'  ENTER  VALUE  OF  HIGH  %  C02 :'//) 
READ(*,61)  CONC2 
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W  =  2  *  CONC2  *  MEANHIGH 
X  =  CONC2  •  (MEAN ZERO  +  MEANHIGH) 

Y  =  2  *  CONC2  •  CONC2 
Z  =  CONC2  *  CONC2 
REGCOEFF  =  (W  -  X)  /  (Y  -  Z) 

V  =  MEANZERO  +  MEANHIGH 
CONSTANT  =  (V  -  REGCOEFF  *  CONC2)  /  2 

c 

WRTTE(*,7)  REGCOEFF, CONSTANT 
WRITE(*,82) 

82  FORMAT(///'  ENTER  NAME  OF  CALIBRATION  FILE:'//) 

WRITE(*,4000)  BELL 
c 

READ(*,73)  FILE2 

OPEN(3,FILE = FILE2.STATUS  =  NEW, FORM  =  FORMATTED) 

WRITE(3,8)  REGCOEFF.CONSTANT 
CLOSE(3) 
c 

WRITE(*,85) 

85  FORMAT(///’  WHAT  DO  YOU  WISH  TO  DO  NOW?'/ 

•  '  1  RETURN  =  CALIBRATE  PNEUMOTACHOGRAPH,'/ 

•  ’  6  RETURN  =  REPEAT  BOTH  LAST  CALIBRATIONS,’/ 

•  '  9  RETURN  =  END  Program.'//) 

READ(*,20)  F 

IF  (F.EQ.6)  GOTO  84 
IF  (F.EQ.1)  CONTINUE 
IF  (F.EQ.9)  GOTO  9 

cccccccccrfffccccfffgccfccccccccccccccccccccccccccccaxcccccccccccccccccccccaxccccccccccccccccccc 
c  Pneumotachograph  calibration 
c 

c  Initialising  arrays  LITER(J)  &  NOFLOW(J). 
c 

86  DO  8600  J  =  1,2002 
LITER(J)  =  0 

8600  CONTINUE 
SAVE=0 


DO  8601  J  =  1,102 
NOFLOW(J)=0 
8601  CONTINUE 
c 

IMODE  =  0 
PARAM(l)  =  #300 
PARAM(2)  =  2 
PARAM(3)  =  3 

RCODE  =  DASH  16(IMODE,  PARAM) 
c 

IMODE  =1 
PARAM(l)  =  4 
PARAM(2)  *  4 
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RCODE  -  DASH16(IMODE,  PARAM) 
c 

IMODE  =  17 
PARAM(l)  =  1000 
PARAM(2)  =  20 

RCODE  =  DASH  16( IMODE,  PARAM) 
c 

WRITE(*,860) 

860  FORMAT(//'  CONNECT  PNEUMOTACHOGRAPH  TRANSDUCER  TO  CHANNEL  FOUR.  ) 
WRITE(*,861) 

861  PORMAT(///'  *•*  PRESS  TRIGGER  TO  COLLECT  BASEUNE  DATA  •••'//) 
c 

IMODE  =19 
PARAM(1)  =  2 
PARAM(2) = 300 
PARAM(3)=0 

RCODE  =  DASH16(IMODE,  PARAM) 
c 

WRITEC31) 

WRfTE(*,4000)  BELL 
c 

IMODE  =  4 
PARAM(l)  =  102 

PARAM(2)=OFFADR(NOFLOW(0)) 

PARAM(3)  =  1 

RCODE  =  DASH16(IMODE,  PARAM) 
c 

WRITEC.4000)  BELL 
4000  FORMATC  >1) 

WRITE(*,4) 

c 

A  =  0 

DO  862  J  =  2,101 
A  =  A  +  NOFLOW(J) 

862  CONTINUE 
c 

c  Writing  to  the  screen  all  the  basal  flow  data,  plus  their  average, 
c 

MEANNOFLOW  =  A  /  100 

WRITE(*,863)  (NOFLOW(J)4  =  2,101), MEANNOFLOW 

863  FORMATflOflOflXJS)/),'  MEAN  BASELINE  FLOW  =  ,15///) 
c 

WRITE(*,42) 

WRITE(*,43) 

READC^O)  G 
IF  (GEQ.l)  CONTINUE 
IF  (G.EQ.6)  GOTO  86 
IF  (G.EQ.9)  GOTO  9 
c 

c  Saving  the  zero  flow  mean  in  a  specific  file, 
c 

WRITEC32) 
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WRITE(*,4000)  BELL 
READ(*,73)  FTLE3 

OPEN(3,FILE  =  FILE3, STATUS  =  'NEW’, FORM  =  FORMATTED) 

WRITE(3,865)  MEANNOFLOW 

865  FORMAT(I7) 

CLOSE(3) 

c 

WRITE(*,866) 

866  FORMAT(///'  CALIBRATE  PNEUMOTACHOGRAPH  USING  1  LITER  SYRINGE.'/ 

*  '  START  FLOW  AS  SUCTION.'/ 

*  '  USE  SIX  (6)  PUMPS  OF  VARIABLE  SPEED.'/ 

*  '  DO  NOT  HIT  THE  PLUNGER  OF  THE  SYRINGE  HARD  AGAINST/ 

*  '  THE  ENDS  OF  ITS  HOUSING,  AS  THIS  PRODUCES  POOR  DATA.'///) 
c 

8651  IMODE  =  0 

PARAM(l)  =  #300 
PARAM(2)  =  2 
PARAM(3)  =  3 

RCODE  =  DASH  16(IMODE,  PARAM) 
c 

IMODE  =  1 
PARAM(l)  =  4 
PARAM(2)  =  4 

RCODE  =  DASH16(IMODE,  PARAM) 

c 

IMODE  =  17 
PARAM(l)  =  1000 
PARAM(2)  =  20 

RCODE  =  DASH  16(IMODE,  PARAM) 

c 

WRITE(Vl) 

c 

IMODE  =19 
PARAM(1)  =  2 
PARAM(2)  =  300 
PARAM(3)  =  0 

RCODE  =  DASH16(IMODE,  PARAM) 
c 

WRITE(*31) 

WRITE(*,4000)  BELL 

c 

IMODE =4 

PARAM(l)  =  2002 

PARAM(2)  =  OFFADR(LITER(0)) 

PARAM(3)  =  1 

RCODE  =  DASH16(IMODE,  PARAM) 
c 

WRITER, 4000)  BELL 
WRITE(*,4) 

WRITE(*,867) 

867  FORMAT^//1  WHAT  DO  YOU  WISH  TO  DO  NOW?'/ 

*  '  1  RETURN  =  SAVE  DATA,'/ 
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*  •  6  RETURN  -  LOOK  AT  DATA.  /) 

READ(*,20)  H 

IF  (H.EQ.1)  GOTO  870 
IF  (H.EQ.6)  GOTO  868 
c 

c  It  is  advisable  to  look  at  the  data  to  check  two  things: 
c  (1)  that  all  syringe  pumps  are  wholly  present;  and 
c  (2)  that  no  data  points  have  flooded  the  A/D  board.  This  occurs 
c  when  the  flow  signal  produces  an  output  greater  than  +  /-  5  Volts, 
c  The  board  only  converts  data  in  this  range.  All  data  outside  this 
c  range  will  be  given  maximal  binary  equivalents  (i.e.  2048  or  -2048). 
c 

868  WRITEC3680)  (UTER(J)4  =  L2000) 

8680  FORMAT  (10(17)) 

c 

864  WRITE(\869) 

869  FORMAT(///'  WHAT  DO  YOU  WISH  TO  DO  NOW?'/ 

*  '  1  RETURN  =  SAVE  DATA,'/ 

*  '  6  RETURN  =  REPEAT  6  LITER  FLOW  CALIBRATION,'/ 

*  '  9  RETURN  =  END  Program.'//) 

READ(*,20)  H 

IF  (H.EQ.1)  CONTINUE 
IF  (H.EQ.6)  GOTO  8651 
IF  (H.EQ.9)  GOTO  9 
c 

870  WRITE(*,82) 

WRITE(*,4000)  BELL 

READ(*,73)  FILE4 

c 

c  Saving  flow  calibration  data  to  a  specific  file, 
c  These  data  are  integrated  to  derive  the  calibration  value 
c  for  the  number  of  binary  units  equivalent  to  1  litre.  An  example 
c  of  this  is  given  following  this  program.  We  typically  do  this  in 
c  the  analysis  program,  but  it  can  be  inserted  here  and  saved  instead 
c  of  the  full  flow  calibration  data.  We  save  the  full  data  in  case  we 
c  get  some  unexpected  errors  which  were  not  apparent  during  calibration, 
c 

OPEN(3,FILE = FILE4, STATUS  =  NEW, FORM  =  UNFORMATTED) 

WRITE(3)  (UTER(jy  =  1,2000) 

SAVE-1 

CLOSE(3) 

WRITE(*,900) 

900  FORMAT(///'  CALIBRATION  COMPLETED,  AND  ALL  DATA  SAVED’/ 

*  ’  WHAT  DO  YOU  WISH  TO  DO  NOW?'/ 

*  ’  1  RETURN  =  GO  TO  START  OF  Program'/ 

*  '  9  RETURN  -  TERMINATE  Program.  //) 

READ(*,20)  ANS 

IF  (ANS.EQ.1)  GOTO  100 
IF  (ANS.EQ.9)  CONTINUE 
c 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
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r 


i 


c  Subroutine  that  prevents  program  termination  without  first 
c  saving  the  flow  calibration  data  to  a  specific  file, 
c 


9 

9989 


IF  (SAVE.LT.1)  GOTO  864 
WRrrEC,9989) 

FORMAT(/ / //// / / 

. 7 

.  . 7 

•  *••**•*•  GOOD  BYE  *••*•*•*•/ 

.  . . 7 

'  . *••*•7///) 

STOP 

END 


cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccco 
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PART  II:  Integration  of  the  calibration  How  signal. 


Integration  subroutine  used  with  the  flow  calibration  data  to  obtain 
the  number  of  computer  /binary  units  that  are  equivalent  to  1  litre  of 
the  volume  calibration  standard.  This  subroutine  may  be  seen  in  the 
composite  program  in  Appendix  7. 


cccccccccccccccccccccrcccccccccfccffcccccccccccccfcccccccccccccccccccccccccccccccccccccccccccccccc 
c  Pre-requisites  for  this  subroutine  include: 
c  (1)  calibration  data  filenames  (NAME1  and  NAME2); 
c  (2)  a  set  flow  limit  -  this  is  the  binary  equivalent  of  typical 
c  flow  encountered  during  the  experiment  (this  is  variable 

c  between  experiments,  and  between  amplifier  and  their  gain 

c  settings); 

c  (3)  a  arbitrary  volume  limit  -  swallowing,  sighing  etc  will  produce 
c  spurious  volume  excursions  -  this  limit  is  used  to  prevent 

c  these  volumes  being  used  durinr  data  analysis,  but  is  also 

c  used  here  to  prevent  artifact  a  om  being  included  in  the 

c  calibration; 

c  (4)  a  volume  for  the  calibration  syringe,  determined  by  repeated 
c  water  filling. 

c 

FILEBASE=NAME1 
FILEFLOV  =  NAME2 
FLOW  =500 
ARB1TLIM  =  500 
SYRINGE  =  3.000 
c 

c  Reading  the  zero  flow  calibration  data  from  data  disk, 
c 

OPEN(3,FILE = FILEB  AS  E,ST ATU  S  =  OLD) 

READ(3,101)  ZEROFLOW 
WRITE(*,102)  ZEROFLOW 

101  FORMAT(I7) 

102  FORMAT('  BASAL  FLOW  (RAW  UNITS)  =  ,15//) 

CLOSE(3) 

50  FORMAT(2XJJ9Z^!X,F92) 

c 

c  Calculation  of  baselimit  used  to  check  for  flow  later 
c  in  program.  This  is  pre-programd  to  suit  anticipated 
c  flow  when  volume  is  sucked  slowly  through  the  pneumotach  to 
c  gve  a  negative  signal.  Transducer  =  Validyne  i  2.0  cmH20. 
c  The  4/100  value  is  used  to  set  a  4%  error  level, 
c  Two  limits  are  set:  (1)  BASELIM  for  inspiration  and  (2)  TBASELIM 
c  for  expiration.  These  serve  as  a  zone  above  and  below  the  actual 
c  zero  flow  level  within  which  the  flow  signal  can  oscillate  before 
c  triggering  integration.  This  is  very  important  in  the  analysis 
c  subroutines  since  it  minimizes  the  effect  of  artifacts  without 
c  compromising  the  integrity  of  the  analysis. 
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c 


BASE  =  (FLOW  -  ZEROFLOW)  *4/100 
BASELIM  =  BASE  +  ZEROFLOW 
TBASELIM  =  ZEROFLOW  -  BASE 
C  WRrTE(*,1345)  BASE 

1345  FORMATO  VALUE  OF  BASE  =  15//) 

WRITE(*,104)  BASELIM 
WRITE(*,104)  TBASELIM 

104  FORMATO  VALUE  OF  BASELIMIT  =  '15//) 
c 

c  Reading  the  flow  calibration  data, 
c 

OPEN(3,FTLE  =  FILEFLOW, STATUS  =  OLD', FORM  =  UNFORMATTED) 
READ(3)  (UTER(jy  =  1,2000) 

11  FORMAT(10(I7)) 

CLOSE(3) 

c 

c  Integration  of  flow  to  derive  volume, 
c  L  will  represent  the  breath  counter  (i.e.  #  of  breaths), 

c  J  will  represent  the  sample  number  currently  analysed, 

c 

c  Checking  whether  flow  is  baseline  or  a  calibration  input, 
c 

M  =  1 
L=0 

12  DO  120  J  =  M,2000 

IF  (LITER(J).LT.BASELIM)  GOTO  1200 

120  CONTINUE 
1200  CONTINUE 
c 

IF  (J.GE.2000)  GOTO  129 
c 

c  Integration  of  calibration  flow  data, 
c 

VOLINSP  =  0.0 
c 

121  DO  127  T  =  J,2000 

I  VOL = UTER(T)-ZEROFLO  W 
R  VOL = FLOAT(IVOL) 

VOUNSP  =  VOUNSP-R  VOL 
IF  (LITER(T).GE.BASELIM)  GOTO  1270 

127  CONTINUE 

1270  CONTINUE 

c 

c  Checking  data  for  false  volumes  before  storing  volumes, 
c 

IF  (VOLINSP.LE-ARBITLIM)  GOTO  128 
L=L+ 1 

CALVOL(L)  =  VOUNSP 

128  IF  (T.GE.2000)  GOTO  129 

M-T 
GOTO  12 
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c 

c  Calculation  of  average  calibration  value  (binary  units  per  litre), 
c 

129  RVOLSUM =0.0 
DO  1280  W=1JL 

RVOLSUM  =  RVOLSUM  +  CALVOL(W) 

1280  CONTINUE 
c 

VOLAVE  =  RVOLSUM  /  (SYRINGE*L) 

1281  FORMAT(10(I7)) 

WRITE(*,1282)  VOLAVE 

1282  FORMATC  MEAN  FLOW  CALIBRATION  VALUE  =\F14.4,'  units/litre'//) 
c 


APPENDIX  3 


Respiratory  pattern  subroutines. 


Respiratory  pattern  subroutines:  These  routines  perform  two  functions: 

(1)  compute  tidal  volumes  from  data  obtained  during  rest  or  exercise; 

(2)  compute  all  respiratory  timing  (including  pauses)  data  and  derive 
minute  ventilation.  These  routines  are  contained  within  the  composite 
example  in  Appendix  7. 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
c  Initializing  variables 
c 

DO  3002  J=  1,40 
TTDALV(J) = 0.0 
VOLBTPS(J) =0.0 
IONSET(J)=0 
EONSET(J)=0 
TINSP(J)  =  0.0 
TEXP(J)=0.0 
TTOTAL(J) =0.0 
3002  CONTINUE 
c 

TOTAL  VOL  =  0.0 

c 

c  COUNTERS:  I  =  inspiration  onset,  E  =  expiration  onset,  [onset 
c  occurs  when  the  flow  signal  crosses  a  limit  threshold] 

c  P  =  pairs  of  insp.  and  exp.,  M  =  temporary  marker, 

c  L  =  tidal  volumes. 

c  NOTE:  I  is  used  as  a  respiratory  cycle  counter  from  this  point, 
c  Integrating  flow  to  obtain  volume  (inspiration  -ive). 
c 

1  =  0 
E=0 
M  =  1 
L=0 
c 

c  Here  we  are  checking  to  make  sure  we  start  the  analysis  with 
c  a  signal  in  the  inspiratory  phase  [convention:  inspiration 
c  produces  a  negative  flow  signal].  The  program  doesn't  move  forward 
c  till  the  desired  point  is  achieved  in  the  flow  signal, 
c 

3  DO  3111  J  =  M,MAX 

IF  (PNEUMO(J).GE.BASELIM)  GOTO  3112 

3111  CONTINUE 

3112  DO  31  K=J,MAX 

IF  (PNEUMO(K).LT.BASELIM)  GOTO  310 
31  CONTINUE 

c 

310  IF  (K.GE.MAX)  GOTO  3102 

c 
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C  WRTTE(*,5500)  K,PNEUMO(K) 

C5500  FORMAT('  IONSET  AT  POINT  '45,'  =  ’,15) 

C 

VOLINSP  =  0.0 
VOL2INSP=0.0 
EXVOLINSP = 0.0 
EX2VOLINSP  =  0.0 
c 

c  Since  we  know  the  signal  is  negative  (above),  we  can  now  look  for 
c  the  next  onset  of  inspiration.  This  ensures  that  the  first  inspiration 
c  is  trully  a  complete  breath  and  not  partially  cut  by  failure  of  the 
c  subject  to  coincide  inspiration  onset  with  the  start  of  data  collection, 
c  Integration  of  inspiratory  flow  (addition,  ie  VOLINSP-(-RVOL)). 
c 

DO  32  T = K,MAX 

IF  (PNEUMO(T).GE.BASELIM)  GOTO  33 
I  VOL = PNEUMO(T)-BASEFLOW 
RVOL=FLOAT(IVOL) 

VOLINSP  =  VOLINSP  -  RVOL 

32  CONTINUE 

33  CONTINUE 
c 

c  The  -RVOL  converts  volume  to  positive, 
c  Check  for  small,  false  volumes, 
c 

c  WRITE(*,5501)  VOLINSP 

C5501  FORMATO  INTEGRATED  TOTAL  =  \F10.2) 

c 

IF  (VOLINSP.LT .ARBITLIM)  GOTO  3103 
c 

c  Marking  points  of  inspiration  start  (IONSET),  &  end  (EONSET), 
c  and  expiration  start  (STARTEXP),  &  end  (ENDEXP). 
c 

1  =  1  +  1 

IONSET(I)  =  K 
EONSET(I)  =  T 
DO  330  J  =  1,15 
VOL(J)  =  0.0 
330  CONTINUE 
c 

L  =  L+ 1 

c 

c  Applying  the  calibration  value  (VOLAVE:  see  Appendix  2). 
c 

VOL(L)  =  VOLINSP  /  VOLAVE 


c 

c  Volume  correction  for  regulator  pressure  and  PLC,  then  BTPS 
c  correction  (all  values  are  in  cmH20  and  degC  units), 
c  Variables:  PREG  =  regulator  pressure  at  the  center  of  the  diaphragm 
c  (if  needed),  PH20  =  water  vapour  pressure  at  TA,  TA  =  ambient  temperature, 
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c  PL C  =  Jung  centroid  pressure  (if  needed),  64.056  =  47.1  *136  (cmH20), 

c  TA2  =  temperature  of  expired  gas  (heater  element  involved),  PH202  = 

c  expired  water  vapour  pressure  at  temp  TA2. 
c 

VOL(L)  =  VOL(L)  *  (PBAR+PREG)  /  (PBAR) 

VOLBTPS(L)  =  VOL(L)*310/(273 + TA)*(PBAR  +  PLC-PH20) 

*  /(PBAR  +  PLC-64.056) 

TIDALV(L)  =  VOLBTPS(L) 

c 

3103  IF  (T.GE.MAX)  GOTO  3102 

IF  (VOLINSP.GTARBITLIM)  GOTO  3104 
M  =  T 
GOTO  3 
c 

c  This  process  is  now  duplicated  for  the  expiratory  phase.  This  phase 
c  may  be  omitted  if  only  inspiratory  data  are  collected, 
c 

3104  DO  34  EX  =  T,MAX 

IF  (PNEUMO(EX)  .GT.TBASELIM)  GOTO  340 
34  CONTINUE 

340  CONTINUE 
c 

DO  341  EXP = EX,  MAX 

IF  (PNEUMO(EXP).LE.TBASELIM)  GOTO  3410 
E  VOL = PNEUMO(EXP)-B  ASEFLOW 
ER  VOL = FLOAT(E  VOL) 

EXVOUSNP  =  EXVOLINSP  +  ERVOL 

341  CONTINUE 

3410  CONTINUE 
c 

IF  (EXP.GE.MAX1  GOTO  3102 
IF  (EXVOLINSP.GT .ARBITLIM)  GOTO  3411 
T  =  EXP 
GOTO  3104 

C 

3411  E  =  E+ 1 
STARTEXP(E)  =  EX 
ENDEXP(E)  =  EXP 

c 

EXTIDALV(L)  =  EXVOLINSP/VOLAVE*(PBAR  +  PREG)/(PBAR)  *310/(273 + TA2) 

*  *(PBAR  +  PLC-PH202)/(PBAR  +  PLC-64.056) 
c 


M  =  T 
GOTO  3 

3102  CONTINUE 
c 

cccccccccccccccccccccccccccccccccccccccccccccccccccceccccccccccccccccecccccccccccccccccccccccccccc 

c 

c  RESPIRATORY  TIMING  SUBROUTINES  (including  post-insp 
c  [POSTINSP]  and  post-expiratory  pauses  [POSTEXP]). 
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c 

cccccccccccceccccccccccccccccccecccecccceccccccccceccecccccceccecccccccecceeccccccccccccccccccc^ 

c 

c  Computing  timing/pausing  periods  based  upon  a  sampling  rate  of 
c  50  Hz.  The  error  is  this  calculation  is  20  milliseconds,  and  can 
c  be  reduced  by  increasing  the  sampling  rate. 

c  Variables:  IONSET  =  start  of  inspiration,  EONSET  =  end  of  inspiration, 
c  STARTEXP  =  start  of  expiration,  END  EXP  =  end  of  expiration,  POSTTNSP  = 
c  post-inspiratory  pause,  POSTEXP  =  post-expiratory  pause 
c 

DO  42  P=1,(I-1) 

TTNSP(P)  =  (EONSET(P)  -  IONSET(P))  /  50.0 
TEXP(P)  =  (ENDEXP(P)-STARTEXP(P))/50.0 
TTOTAL(P)  =  (IONSET(P  +  1)-IONSET(P))/50.0 
POSTINSP(P)  =  (STARTEXP(P)-EONSET(P)) /50.0 
POSTEXP(P)  =  (IONSET(P+1)-ENDEXP(P))/50.0 
42  CONTINUE 

c 

WRTTE(*,103)  I,E 

103  FORMAT('  ,15, 1  INSPIRATIONS  AND  ,13,'  EXPIRATIONS'///) 

WRITE(*,420) 

420  FORMAT('  TIDAL/I  TIDAL/E  INSP/Time  EXP /Time  T/tot ' 

*  '  POST-INSP/Pause  POST-EXP /Pause) 

c 

DO  421  P  =  1,(I-1) 

WRITE(*,422)  TIDAL V(P),EXTIDALV(P),TINSP(P),TEXP(P),TTOTAL<P), 

*  POSTINSP(P),POSTEXP(P) 

421  CONTINUE 

422  FORMAT(F7.3,2X,F7.3,2X,F7.4,lX,F7.4,lX,F7.4,5X,F7.4,8X,F7.4) 
c 

AVTINSP =0.0 
AVTEXP=0.0 
A  VTTOTAL =0.0 
AVTIDALV  =  0.0 
AVINSPPOST  =  0.0 
AVEXPPOST  =  0.0 
c 

c  Summing  the  timing  components 
c 


DO  44  P=  1,(I-1) 

AVTIDALV  =  AVTIDALV  +  TIDAL  V(P) 
AVTINSP  =  AVTINSP  +  TINSP(P) 

AVTEXP  =  AVTEXP  +  TEXP(P) 

A  VTTOTAL  =  A  VTTOTAL  +  TTOTAL(P) 
AVINSPPOST  =  AVINSPPOST  +  POSTINSP(P) 
AVEXPPOST  =  AVEXPPOST  +  POSTEXP(P) 
44  CONTINUE 
c 

c  Computing  averages  for  each  timing  component. 
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AVTIDALV  =  AVTIDALV  /  (1-1) 

AVTINSP  =  AVTINSP  /  (1-1) 

AVTEXP  =  AVTEXP  /  (1-1) 

AVTTOTAL  =  AVTTOTAL  /  (1-1) 

RFREQ = 60.0/A  VTTOTAL 
VTTIMEINSP  =  AVTIDALV  /  AVTINSP 
MINVENT  =  AVTIDALV  *  RFREQ 
TINTOTAL  =  AVTINSP  /  AVTTOTAL 
AVINSPPOST  =  AVINSPPOST  /  (1-1) 

AVEXPPOST  =  AVEXPPOST  /  (1-1) 

c 

WRITE(*,440)  AVTIDALV 

440  FORMAT(//////'  AVERAGE  TIDAL  VOLUME  =  '.F12.4/) 

WRITE(*,441)  AVTINSP 

441  FORMAT('  AVERAGE  TIME  INSP  =  ,F12.4) 

WRITE(*,442)  AVTEXP 

442  FORMATC  AVERAGE  TIME  EXPIR.  =  \F12.4) 

WRITE(*,443)  AVTTOTAL 

443  FORMAT('  AVERAGE  TOTAL  TIME  =  ,F12.4/) 

WRITE(*,444)  VTTIMEINSP 

444  FORMATC  VT/Ti  =  .F12.4) 

WRITE(*,446)  TINTOTAL 
WRITE(  *  ,4400)  RFREQ 

4400  FORMATC  BREATHING  FREQUENCY  =  \F12.4) 

WRITE(*,445)  MINVENT 

445  FORMATC  MINUTE  VENT.  =  ,F12.4/) 

446  FORMATC  Ti/Ttot  =  '^12.4/) 

WRITEC,447)  AVINSPPOST 

447  FORMATC  AVERAGE  POST-INSPIRATORY  PAUSE  =  \F12.4) 

WRITE(*,448)  AVEXPPOST 

448  FORMATC  AVERAGE  POST-EXPIRATORY  PAUSE  =  \F12.4) 

c 

CCCCCCCCCCCCCCCCCCCCOXCCCC£CCCCOQCC*XCCCCXX£CCCCCCCCCC£UXCCCCCCCCCaX£CC£CCCCCCCCCCCCCCCCCCCCCCC 
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APPENDIX  4 


The  determination  of  the  various  lung  volumes  is  first  dependent  upon  the 
analyses  performed  in  Appendix  3,  since  each  of  the  volumes  and  capacities 
shall  be  considered  as  a  separate  tidal  volume  (TIDAL V).  This  subroutine  is  contained 
within  the  example  program  in  Appendix  7. 


c 

WRITEC3337) 

DO  3339  P=  LI 

WRITEC^Sg)  P,  TTDALV(P) 

3339  CONTINUE 

3337  FORMATO  VOLUME  #  SPIROMETRIC  VOLUMES  (inspiratory)') 

3338  FORMAT(2X,I4,UX,F73) 


APPENDIX  5 


Maximal  voluntary  ventilation  subroutine 


Calculations  performed  here  are  also  dependent  upon  those  performed  in  Appendix  3,  since  each  tidal 
volume  (TIDAL V)  is  added  for  the  duration  of  the  maneuver,  and  the  duration  from  the  start  to  the  end  of 
the  maneuver  are  required  (MWTIME).  This  subroutine  is  contained  within  the  example  program  in 
Appendix  7. 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 

MW  =  0.0 
c 

DO  35  P=  1,(I-1) 

MW  =  MW  +  TIDAL  V(P) 

35  CONTINUE 

c 

c  The  number  50.0  and  60.0  are  used  because  the  sampling  rate  was  50  Hz  and  the 
c  output  is  required  in  liters  per  minute  (60.0). 
c 

MWTIME  =  (ENDEXP(I-l)-IONSET(l))  /  50.0  /  60.0 
MW  =  MW  /  MWTIME 
WRITE(*,350)  MW 

350  FORMAT(///'  MAXIMAL  VOLUNTARY  VENTILATION  FOR  THIS'/ 

*  '  TRIAL  =  JF12.4/  Litres  BTPS'///) 

c 

cccccccccccccccccccccccccaxcccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
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APPENDIX  6 


Maximal  flow-volume  loop  subroutines. 


Calculations  performed  here  are  also  dependent  upon  those  performed  in  Appendix  3,  since  the  onset  of 
each  of  the  phases  of  respirations  are  used  as  markers  for  the  analysis.  This  subroutine  is  contained  within 
the  example  program  in  Appendix  7.  Several  variables  have  previously  been  explained  in  other  Appendices 
and  their  explanations  are  not  repeated. 


cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 

c  Convention  for  flow  signs  during  analysis:  inspiration  -ive, 
c  expiration  +ive. 
c 

c  In  these  routines  all  respiratory  cycles  are  analysed,  thus 
c  in  DO  6  P-1,I  and  not  1-1.  Requires  a  data  output  filename  (FILE2). 
c  TOPVOL  =  the  final  lung  volume  or  TLC,  RVOLIN  =  integrated  inspired 
c  volume,  VOLEX  =  integrated  expired  volume,  RINFLOW  =  inspiratory  flow 
c  variable,  R EXFLOW  =  expiratory  flow  variable, 
c 

OPEN(4,FILE = FILE2,STATUS  =  NEW1) 

TOPVOL =0.0 
c 

DO  6  P  =  1,1 
c 

c  Setting  arrays  and  variables  to  zero, 
c 

DO  60  J=  1,200 
RVOLIN(J)  =0.0 
VOLEX(J)=0.0 
RINFLOW(J)=0.0 
R  VOLEX(J) = 0.0 
REXFLOW(J)=0.0 
60  CONTINUE 

c 

c  Computes  the  inspiratory  volume  and  flow  at  each  20  millisecond 
c  interval.  These  data  are  written  to  both  the  screen  and  output  file, 
c 

WRITE(*,63) 

DO  610  PPP  =  IONSET(P),(EONSET(P)- 1) 

PP  =  (PPP-IONSET  (P))  + 1 

INVOL = PNEUM  O(PPP)  -  B  AS  EFLO  W 

RINVOL=  FLOAT(INVOL) 

RINVOL=RINVOL/VOLAVE*(PBAR  +  PREG)/(PBAR)*310/ 

*  (273  +  TA)*(PBAR  +  PLC-PH20)/(PBAR  +  PLC-64.056) 

IF  (RINVOL.LT. (0.0))  GOTO  6100 
RINVOL = 0.0-RINVOL 


41 


6100  RVOLIN(PP)«RVOLIN(PP-l)-RINVOL 

RINFLO  W(PP)  =  FLOAT  (INVOL) 

RINFLOW(PP)  =  RINFLOW(PP)/VOLAVE*50.0*(PBAR  +  PREG)/(PBAR) 

*  *310/(273 + TA)*(PBAR  +  PLC-PH20)/(PBAR  +  PLC-64.056) 

WRITE(*,631)  PPP,RVOLIN(PP),RINFLOW(PP) 

WRITE(4,631)  PPP,RVOUN(PP),RINFLOW(PP) 

610  CONTINUE 

TOPVOL = RVOLIN(PP) 

WRITE(*,630) 

c 

c  Expiratory  data  analysis, 
c 

DO  620  PPP  =  STARTEXP(P),(ENDEXP(P)-1) 

PP = (PPP-STARTEXP(P))  + 1 
EVOL = PNEUMO(PPP)-BASEFLOW 
REXVOL = FLOAT(EXVOL) 

REXVOL = REXVOL/VOLAVE*(PBAR  +  PREG)/(PBAR)*310/ 

*  (273+TA2)*(PBAR  +  PLC-PH202)/(PBAR  +  PLC-64.056) 

RVOLEX(PP)  =  RVOLEX(PP-l)  +  REXVOL 

VOLEX(PP) = TOPVOL-R  VOLEX(PP) 

REXFLOW(PP)  =  FLOAT(EXVOL) 

REXFLOW(PP)  =  REXFLOW(PP)/VOLAVE*50.0*(PBAR  +  PREG)/(PBAR) 

*  *310/(273 +TA2)*(PBAR  +  PLC-PH202)/(PBAR  +  PLC-64.056) 

WRITE(*,631)  PPP, VOLEX(PP) .REXFLOW (PP) 

WRITE(4,631)  PPP,VOLEX(PP),REXFLOW(PP) 

620  CONTINUE 
6  CONTINUE 
c 

63  FORMATC  INSPIRATORY  DATA'/ 

*  '  #  Vol(BTPS)  Row  ') 

630  FORMAT(/'  EXPIRATORY  DATA'/ 

*  '  #  Vol(BTPS)  Row  ') 

631  FORMAT(  150(I4,1X,F8.4,2X,F9.4) ) 

CLOSE(4) 

c 

c  Subroutines  for  computing  FEV1.0  from  the  data  obtained 
c  during  the  maximal  flow-volume  loop  trials, 
c 

EVOL=0 
FEV1-0.0 
EXVOLINSP =0.0 

DO  64  P  =  STARTEXP(1),(STARTEXP(  1)  +  49) 

EVOL=  PNEUMO(P)-BASEFLOW 
ER  VOL=  FLOAT(EVOL) 

EXVOLINSP  =  EXVOLINSP  +  ERVOL 

64  CONTINUE 
c 
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FEV1  -  EXVOUNSP/VOLAVE*(PBAR  +  PREG)/(PBAR)*310 
*  /(273 + TA2)*(PBAR  +  PLC-PH202)/(PBAR  +  PLC-64.056) 
WRITE(*,65)  FEV1 

65  FORMAT(///'  FEV1  =  ',F8.4/  Litres  BTPS.’///) 
c 


APPENDIX  7 
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Many  of  the  comments  necessary  to  understand  the  following  program  and  its  variables  have  been  omitted, 
since  they  have  been  incorporated  in  the  Appendices  above. 


c 

c  PROGRAM  TO  ANALYSE  DATA  FROM  RESPIRATORY  VOLUME  EXPERIMENTS, 
c 

c  Program  uses  data  obtained  from  A/D  program  SPIRO-D.EXE 
c  using  the  DASH16  A/D  board  (MetraByte);  (50  Hz,  2  channels) 
c 

c  INPUT  DATA  REQUIRED:  Ch2  =  Trigger,  Ch3  =  flow 

c  Flow  signal  is  both  inspiratory  and  expiratory 
c 

c  NOTE:  THIS  Program  HAS  OUTPUT  IN  S.I.  UNITS 

c  THE  ANALYSIS  RELIES  UPON  USING  A  HEATED  PNEUMO  @  37oC 
c 

c  THE  FOLLOWING  ANALYSES  ARE  PERFORMED: 

c  (1)  Tidal  volumes,  (2)  Respiratory  timing  (Ti,  Te,  Ttot, 
c  Vt/Ti,  Ti/Ttot,  fb,  Ve,  and  pauses),  (3)  Maximal  Voluntary 
c  Ventilation,  (4)  Flow- Volume  loops  and  FEV1,  (5)  Spirometry 
c  for  determining  lung  volumes  and  capacities, 
c 

c  NIGEL  A.S.  TAYLOR.  MAY  1989 

c  Diving  Life  Support  Equipment  Program 

c  Naval  Medical  Research  Institute 

c  Bethesda,  MD. 

c 

c  LATEST  UPDATE:  4/23/90 
c 

c  This  program  allows  for  unequal  inspired/expired  volumes 
c  by  using  only  complete  cycles,  starting  with  an  inspiration, 

c  Inspiration  is  detected  from  a  negative  signal  from  the  flow 

c  channel, 
c 

c  Program  is  able  to  analyse  data  obtained  when  subjects  are  immersed 
c  and  air  is  pro\ .ed  at  various  hydrostatic  pressures.  This  requires 
c  modification  after  completion  of  the  calibration  analysis  routines, 
c 


c 

$STORAGE:2 

c 

INTEGER  ZEROFLOW, LITER, BASELIM,PHASEl,PHASE2,INVOL 
INTEGER  PNEUMO,NOFLOW,FLOW,BASEIDATA,COMP,EXVOL 
INTEGER  GAIN,PREG,PSTERN,IONSET,EONSET,TIME>lAX 
INTEGER  MAX2,IVOL,BASEFLOW,PRESS, INFLOW, EXFLOW 
INTEGER  A3,C^,E/,U4U^^,PWPP,Q,R&T,W,X,Y,Z 


c 


INTEGER  TBASELIM,EX,EXP, END  EXP, STAR  TEXP,EVOL,REPEAT 


REAL  VOLUME, V0LAVE,RV0LSUM,PH20,TAtAVTIDALV,TINSP 
REAL  TOTALVOL,INSPVOL,SYRINGE,VOL,TIDALV,TEXP,TTOTAL 
REAL  VOLINSP,CALVOL,RVOLARBITLIM,VTTIMEINSPAVnNSP 
REAL  AVTEXPAVTTOTAL,RFREQ,MINVENT,MINV2,TINTOTALTIDE 
REAL  ERVOL,MW,MWTIME,POSTINSP,POSTEXP,TOPVOL,TA2 
REAL  COMPRESS, PLC,SIPLC,REXV0L,V0LEX,RV0LEX,PH202 
REAL  RINVOL,RVOLIN,RINFLOW,REXFLOW,VOL2INSP,EX2VOLINSP 
REAL  INT^XTEXnDALV^XSMALLTV.EXTIDE 
REAL  EXVOLINSP,AVINSPPOST,AVEXPPOST,FEVl 
REALM  PBAR 

DIMENSION  ID  AT A(2,6000)  4-ITER  (2000) 

DIMENSION  PNEUMO(30(X)),VOL(40),VOLBTPS(40),IONSET(40),TINSP(40) 
DIMENSION  CALVOL(10),EONSET(40),TIDALV(40),TEXP(40),TTOTAL{40) 
DIMENSION  POSTINSP(40),POSTEXP(40),PINSP(40) 

DIMENSION  RVOLiN(200),RINFLOW(200),RVOLEX(200),VOLEX(200) 
DIMENSION  REXFLOW(200),TIDE(40) 

DIMENSION  EXTIDALV(40),EXSMALLTV(40),EXTIDE(40) 

DIMENSION  ENDEXP(40),STARTEXP(40) 

CHARACTER  BELL 

CHARACTERS  FILEBASE, FILE  FLO  W,FILEI,FILE2 
CHARACTER  *40  NAME1,NAME2,NAME3,NAME4 
BELL  =  CHAR(7) 

OPENC) 


c 

:  INTEGRATION  OF  FLOW  TO  GET  TIDAL  VOLUMES. 

c 


REPEAT  =  0 
TIME  =  0 
WRITE(*,7000) 

7000  FORMAT(//////// 

•  '  RESPIRATORY  VOLUME  ANALYSIS  Program:'//// 

•  '  DO  YOU  WISH  TO  USE  DEFAULT  CALIBRATION  FILES?'/ 

•  '  1  RETURN  =  YES,'/ 

•  '  9  RETURN  =  NO.'/) 

READ(*,100)  A 

100  FORMAT(I3) 

IF  (A.NE.l)  GOTO  7001 
NAME1  =  VOLBASE.DAT 
NAME2  « •VOLCAUB.DAT 
NAME3  -  VOL2BASE.DAT 
NAME4  =  VOL2CAL.DAT 
GOTO  7002 

7001  WRITE(*,7003) 

7003  FORMAT(//'  ENTER  FOUR  (4)  CALIBRATION  FILENAMES.'///) 


WRITE(  *  ,7004) 

7004  FORMATC  ENTER  LOW  GAIN  BASAL  CALIBRATION  FILENAME:'/) 
READ(\7005)  NAME1 

7005  FORMAT(A40) 

WRITE(*,7006) 

7006  FORMATC  ENTER  LOW  GAIN  FLOW  CALIBRATION  FILENAME:'/) 
READ(*,7005)  NAME2 

WRITE(*,7007) 

7007  FORMATC  ENTER  HIGH  GAIN  BASAL  CALIBRATION  FILEMANE:’/) 
READ(*,7005)  NAME3 

WRITE(*,7008) 

7008  FORMATC  ENTER  HIGH  GAIN  FLOW  CALIBRATION  FILENAME:'/) 
READ(*,7005)  NAME4 

7002  CONTINUE 
1310  WRITE(M) 

1  FORMAT(/// 

*  ’  WHICH  TYPE  OF  DATA  DO  YOU  WISH  TO  ANALYSE?  / 

*  ’  1  RETURN  =  SPONTANEOUS  RESPIRATION/ 

*  '  3  RETURN  =  LUNG  VOLUME  SPIROMETRY/ 

*  '  6  RETURN  =  MAXIMAL  VOLUNTARY  VENTILATION-/ 

*  '  9  RETURN  =  MAXIMAL  FLOW- VOLUME  LOOPS'/) 

READ(*,100)  A 

C 

c  Determining  data  file  dimensions  and  setting  basal  input 
c  according  to  the  gain  settings  used  for  different  experiments 
C 

IF  (A.NE.1AND7V.NEJ)  GOTO  1991 
IF  (A.EQ.1)  TIME =60*50*2 
IF  (A.EQ3)  TIME =45*50*2 
FILEBASE  =  NAME1 
FILEFLOW = NAME2 
FLOW =-300 
ARBITLIM  =  1500 
GOTO  1990 

1991  IF  (A.EQ.6)  TIME  =15*50*2 
IF  (A.EQ.9)  TIME  =12*50*2 
FILEBASE  =  NAM E3 
FILEFLOW  =  NAME4 
FLOW =-500 
ARBITUM  =  500 
1990  MAX2=TIME 
MAX  =  MAX2/2 
C 

SYRINGE  =  3.000 

OPEN(3,FILE  =  FILEBASE,STATUS  =  OLD) 

READ(3,101)  ZEROFLOW 
WRITE(*,102)  ZEROFLOW 

101  FORMAT(I7) 

102  FORMATC  BASAL  FLOW  (RAW  UNITS)  =  ,15//) 

CLOSE(3) 

50  FORMAT(2X,F9.2v2X,F9.2) 
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BASE  =  (FLOW  -  ZEROFLOW)  *4/100 
BASELIM  =  BASE  +  ZEROFLOW 
TBASELIM  -  ZEROFLOW  -  BASE 
WRTTE(*,1345)  BASE 

1345  FORMAT(’  VALUE  OF  BASE  =  T5 //) 

WRITE(*,104)  BASELIM 
WRITE(*,104)  TBASELIM 

104  FORMATC  VALUE  OF  BASELIMIT  =  15//) 
c 

OPEN(3,FILE = FILEFLOW, STATUS  =  OLD', FORM  =  UNFORMATTED) 

READ(3)  (LITER(jy  =  1,2000) 

11  FORMAT(10(I7)) 

CLOSE(3) 

c 

M  =  1 
L=0 

12  DO  120  J  =  M, 2000 

IF  (LITER(J).LT.BASELIM)  GOTO  1200 

120  CONTINUE 
1200  CONTINUE 
c 

IF  (J.GE.2000)  GOTO  129 
VOLINSP=0.0 
c 

121  DO  127  T =J,2000 

IVOL=  UTER(T)-ZEROFLOW 
R  VOL = FLOAT(FVOL) 

VOLINSP=  VOUNSP-RVOL 
IF  (LITER(T).GE.BASELIM)  GOTO  1270 

127  CONTINUE 
1270  CONTINUE 
c 

IF  (VOLINSP.LEARBITLIM)  GOTO  128 
L=L+ 1 

CALVOL(L)  =  VOLINSP 

128  IF  (T.GE 2000)  GOTO  129 
M=T 

GOTO  12 
c 

129  RVOLSUM=0.0 
DO  1280  W=1,L 

RVOLSUM  =  R  VOLSUM  +  CALVOL(W) 

1280  CONTINUE 
c 

VOLAVE  =  RVOLSUM  /  (SYRINGE’L) 

1281  FORMAT(10(I7)) 

WRITE(*,1282)  VOLAVE 

1282  FORMATC  MEAN  FLOW  CALIBRATION  VALUE  =',F14.4,'  uniU/litre’//) 
c 

ccccccccccccc^ccccccce^ccccaxccccccctxcccccccccccccccc^ccccccaxcccacccccccccccccccccccccccca 

c 

IF  (REPEAT.EQ.1)  GOTO  1311 
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c 

2  FORMAT(//'  PLEASE  NAME  THE  DATA  FILE  TO  BE  ANALYSED.’/// 

*  ’  examples:  (1)  BARNES. RES  =  SPONTANEOUS  BREATHING'/ 

*  •  (3)  BARNES.SPI  =  LUNG  VOLUME  SPIROMETRY’/ 

•  •  (6)  BARNES. MW  =  MAXIMAL  VOLUNTARY  VENTILATION1/ 

•  •  (9)  BARNES. FLO  =  MAXIMAL  FLOW- VOLUME  LOOPS’///) 

20  FORMAT(A40) 

c 

WRITE(*,130) 

130  FORMAT('  WHAT  WAS  THE  BAROMETRIC  PRESSURE  (mmHg  F8.4)?'/) 
READC.131)  PBAR 

131  FORMAT(F123) 

PBAR  =  PBAR  *  136 

c 

c  The  question  below  is  used  only  if  inspired  air  temp  was  not  37  °  C 
c  TA  =  inspired  gas  temperature,  TA2  =  expired  temp,  PH20  =  water 
c  vapour  pressure,  PH202  =  expired  water  vapour  pressure.  Both  pressures 
c  are  converted  to  cmH20. 
c 

C  WRITE(*,132) 

C132  FORMATO  WHAT  WAS  THE  INSPIRED  AIR  TEMPERATURE  (deg  F  F8.4)?'/) 

C  READ(*,131)  TA 

TA  =  37.0 
TA2  =  37.0 
PH202  =  47.1  *  136 
WRITE(*,133) 

133  FORMAT(’  WHAT  WAS  INSPIRED  WATER  VAPOR  PRESSURE  (mmHg  F8.4)?'/) 

READ(*,131)  PH20 
PH20  =  PH20  *  136 
c 

c  Variables  used  when  air  is  provided  at  various  breathing  pressures, 
c  COMP  is  compensator  pressure  used  to  balance  the  differential  transducer, 
c  This  is  necessary  when  mouth  pressures  during  immersion  exceed  the  range 
c  limits  of  the  transducer.  PLC  =  lung  centroid  pressure,  COMPRESS  and 
c  SIPLC  are  SI  equivalents  of  the  above  pressures.  PREG  =  hydrostatic 
c  pressure  at  center  of  the  regulator  (ie.  depth  in  cmH20),  PSTERN  = 
c  hydrostatic  pressure  at  the  depth  of  the  sternal  notch,  PHASE  =  the 
c  phase  shift  between  flow  and  pressure  signals  (not  used), 
c 

COMP  =  0 
COMPRESS  =  0.0 
PLC  =  0 
SIPLC  ■=  0 
PREG  =  0 
PSTERN  =  0 
PHASE  =  0 
c 

1311  WRITE(V) 

READ(*,20)  FILE1 
c 

c  The  output  file  is  in  ASCII  format,  and  with  the  extension  .PRN 
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c  it  may  be  IMPORTED  directly  into  a  LOTUS  123  worksheet.  This  is 
c  particularly  useful  when  analysing  the  flow-volume  loops, 
c 

WRITER, 1340) 

1340  FORMAT('  WHAT  DO  YOU  WISH  TO  CALL  THE  OUTPUT  FILE?  // 

•  1  example:  BARN-RES.PRN  ’///) 

READC.20)  FILE2 

c 

135  FORMAT(I3) 
c 

c  Reading  unformatted  data  file 
c 

OPEN(3,FILE  =  FILE1, STATUS  =  OLD, FORM  =  UNFORMATTED) 

DO  21  K=l,2 

READ(3)  (IDATA(K,L),L=  1,MAX2) 

21  CONTINUE 
CLOSE(3) 

c 

c  Placing  flow  data  into  the  array  PNEUMO. 
c 

DO  22  J  =  1,MAX 
PNEUMO(J)  =  IDATA(2rI) 

22  CONTINUE 

24  FORMAT (////'  DATA  FILE  ANALYSED  =  ’A40 ////) 

c 

c  Volume  analysis  routine. 

c  Check  baseline  again,  and  use  same  baselimit  value, 

c  These  next  statements  are  used  to  check  exp.  data  for  basal 

c  points, 
c 

N  =  0 

DO  30  J  =  5,14 
N  =  N  +  PNEUMO(J) 

30  CONTINUE 

NOFLOW  =  N  /  10 
C  WRITEC301)  NOFLOW 

C301  FORMATO  BASAL  FLOW  FROM  TRIAL  (raw  units)  =  ’,15/) 
c 

WRITE(*355)  BELL 
355  FORMATO  ,A1) 

WRITE(*3000) 

3000  FORMAT(///'  Activate  printer  if  print-out  required'/// 

*  ’  1  RETURN  =  TO  CONTINUE/) 

READ(*,100)  F 

WRITE(\24)  FILEl 
c 

BASEFLOW  =  ZEROFLOW 
c 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCXXCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

C 


DO  3002  J  =  1,40 
TIDALV(J)=0.0 
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VOLBTPS(J)  =0.0 
IONSET(J)=0 
EONSET(J)=0 
TINSP(J)=0.0 
TEXP(J)=0.0 
TTOTAL(J) = 0.0 
3002  CONTINUE 
c 

TOTALVOL  =  0.0 
c 

1=0 
E=0 
M  =  1 
L=0 
c 

3  DO  3111  J  =  M,MAX 

IF  (PNEUMO(J)  .GE.BASELIM)  GOTO  3112 

3111  CONTINUE 

3112  DO  31  K=J,MAX 

IF  (PNEUMO(K).LT.BASELIM)  GOTO  310 
31  CONTINUE 

c 

310  IF  (KLGE.MAX)  GOTO  3102 

c 

C  WRITEC4500)  K,PNEUMO(K) 

C5500  FORMAT('  IONSET  AT  POINT  ,15,’  =  ,15) 

C 

VOLINSP=0.0 
VOL2INSP = 0.0 
EXVOLINSP = 0.0 
EX2VOLINSP = 0.0 
c 


DO  32  T  =  K.MAX 

IF  (PNEUMO(T).GE-BASELIM)  GOTO  33 
IVOL=  PNEUMO(T)-BASEFLOW 
RVOL=  FLOAT(IVOL) 

VOLINSP  =  VOLINSP  -  RVOL 

32  CONTINUE 

33  CONTINUE 
c 

C  WRITE(*,5501)  VOLINSP 

C5501  FORMATC  INTEGRATED  TOTAL  =  ’.F10.2) 

c 

IF  (VOLINSP.LT ARBITLIM)  GOTO  3103 
c 

1  =  1  +  1 

IONSET(I)  =  K 
EONSET(I)  =  T 
DO  330  J  =  l,15 
VOL(J)=0.0 
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330  CONTINUE 
c 

L=L+1 

VOL(L)  =  VOLINSP  /  VOLAVE 
c 

c  Volume  correction  for  regulator  pressure  and  PLC,  then  BTPS  correction 
c  (all  values  arc  in  cmH20  and  degC  units), 
c 

VOL(L)  =  VOL(L)  •  (PBAR  +  PREG)  /  (PBAR) 

VOLBTPS(L)  =  VOL(L)*310/(273 + TA)*(PBAR  +  PLC-PH20) 

*  /(PBAR  +  PLC-64.056) 

TIDAL V(L)  =  VOLBTPS(L) 
c 

3103  IF  (T.GE.MAX)  GOTO  3102 

IF  (VOUNSP.GTARBITLIM)  GOTO  3104 
M  =  T 
GOTO  3 
c 

3104  DO  34  EX = T,MAX 

IF  (PNEUMO(EX).GT.TBASELIM)  GOTO  340 
34  CONTINUE 

340  CONTINUE 
c 

DO  341  EXP  =  EX,  MAX 

IF  (PNEUMO(EXP).LE.TBASELIM)  GOTO  3410 
EVOL=  PNEUMO(EXP)-BASEFLOW 
ERVOL= FLOAT(EVOL) 

EXVOUSNP  =  EXVOLINSP  +  ERVOL 

341  CONTINUE 

3410  CONTINUE 
c 

IF  (EXP.GE.MAX)  GOTO  3102 
IF  (EXVOUNSP.GTARBITLIM)  GOTO  3411 
T  =  EXP 
GOTO  3104 
c 

3411  E  =  E+1 
STARTEXP(E)  *  EX 
ENDEXP(E1  =  EXP 

c 

EXTIDALV(L)  =  EXVOUNSP/VOLAVE*(PBAR  +  PREG)/(PBAR)*310/(273+ TA2) 
*  *(PBAR  +  PLC-PH202)/(PBAR  +  PLC-64.056) 

c 

M=T 

GOTO  3 

3102  CONTINUE 
c 

IF  (A.EQ.9)  GOTO  6000 
c 

c  Maximal  voluntary  ventilation  calculations 
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MW  =  0.0 

IF  (A.NE.6)  GOTO  3333 
DO  35  P=1,(I-1) 

MW  =  MW  +  TIDALV(P) 

35  CONTINUE 

MWTIME  =  (ENDEXP(I-l)-IONSET(l))  /  50.0  /  60.0 
MW  =  MW /MWTIME 
WRITE(*r350)  MW 

350  FORMAT(///'  MAXIMAL  VOLUNTARY  VENTILATION  FOR  THIS1/ 

*  •  TRIAL  =  ',F12.4,'  Utres  BTPS '///) 

GOTO  9000 
3333  CONTINUE 
c 

OPEN(4JTLE = FILE2,  STATUS  =  'NEW) 

8000  FORMAT('  ANALYSIS  OUTPUT  FROM  DATA  FILE:  >15/) 

4  FORMAT(*  START /INSP  END/INSP  TIDAL  VOLUME  START/EXP  END/EXP) 

c 

c  Routine  for  writing  spirometric  volumes 
c 

IF  (A.NEJ)  GOTO  3009 
WRITE(*3337) 

DO  3339  P=  1,1 

WRITE(*3338)  P,  TIDAL V(P) 

3339  CONTINUE 

3337  FORMAT('  VOLUME#  SPIROMETRIC  VOLUMES  (inspiratory)) 

3338  FORMAT(2X,I4,UX,F73) 

GOTO  9000 

3009  CONTINUE 
c 

DO  40  P=1,(I-1) 

c  WRITE(4,41)  IONSET(P),EONSET(P),TIDALV(P),STARTEXP(P),ENDEXP(P) 
c  WRITE(*,41)  IONSET(P),EONSET(P),TTDALV(P),STARTEXP(P),ENDEXP(P) 
c40  CONTINUE 

41  FORMAT(15(I6,7X,I63X,F12.4,2X,I6,4X,I6)) 

c 


c  RESPIRATORY  TIMING  SUBROUTINES  (including  post-insp 

c  [POSTINSPj  and  post -expiratory  pauses  [POSTEXP]). 


DO  42  P=  1,(I-1) 

TINSP(P)  -  (EONSET(P)  -  IONSET(P))  /  50.0 
TEXP(P)  =  (ENDEXP(P)-STARTEXP(P))/50.0 
TTOTAL(P)  =  (IONSET(P  +  1)-IONSET(P))/50.0 
POSTINSP(P)  =  (STARTEXP(P)-EONSET(P))/50.0 
POSTEXP(P)  =  (IONSET(P+  l)-ENDEXP(P))/50.0 
42  CONTINUE 
c 

WRITE(M03)  1^ 

103  FORMATC  '45,'  INSPIRATIONS  AND  ',13,'  EXPIRATIONS  ///) 


WRITE(4,420) 

WRITER  *,420) 

420  FORMAT(*  TIDAL/I  TIDAL/E  INSP/Time  EXP /Time  T/tot  ’ 

•  *  POST-INSP /Pause  POST-EXP/Pause  ) 
c 

DO  421  P=1,(I-1) 

WR1TE(*,422)  TIDALV(P),EXTIDALV(P),TINSP(P),TEXP(P),TTOTAL(P), 

•  POSTINSP(P),POSTEXP(P) 

WRrTE(4,422)  TTDALV(P),EXTIDALV(P),TINSP(P),TEXP(P),TTOTAL(P), 

•  POSTINSP(P),POSTEXP(P) 

421  CONTINUE 

422  FORMAT(F7J^X,F73^X,F7.4,lX,F7.4,lX,F7.4^X,F7.4,8X,F7.4) 
c 

AVTTNSP=0.0 
AVTEXP=0.0 
A  VTTOTAL = 0.0 
AVTIDALV  =  0.0 
AVINSPPOST  =  0.0 
AVEXPPOST  =  0.0 
DO  44  P=1,(I-1) 

AVTIDALV  =  AVTIDALV  +  TIDAL  V(P) 

AVTINSP  =  AVTTNSP  +  TINSP(P) 

AVTEXP  -  AVTEXP  +  TEXP(P) 

A  VTTOTAL  =  A  VTTOTAL  +  TTOTAL(P) 

AVINSPPOST  =  AVINSPPOST  +  POSTINSP(P) 

AVEXPPOST  =  AVEXPPOST  +  POSTEXP(P) 

44  CONTINUE 
c 

AVTIDALV  =  AVTIDALV  /  (1-1) 

AVTINSP  =  AVTINSP  /  (1-1) 

AVTEXP  =  AVTEXP  /  (1-1) 

A  VTTOTAL  =  A  VTTOTAL  /  (1-1) 

RFREQ = 60.0/A  VTTOTAL 
VTOMEINSP  =  AVTIDALV  /  AVTINSP 
MINVENT  =  AVTIDALV  *  RFREQ 
TINTOTAL  =  AVTINSP  /  A  VTTOTAL 
AVINSPPOST  =  AVINSPPOST  /  (1-1) 

AVEXPPOST  =  AVEXPPOST  /  (1-1) 
c  MINV2  *  VTTIMEINSP  *  TINTOTAL 
c 

WRITE(*,440)  AVTIDALV 
WRITE(4,440)  AVTIDALV 

440  FORMAT (//////'  AVERAGE  TIDAL  VOLUME  =  \F12.4/) 

WRITE(*,441)  AVTINSP 

WRITE(4,441)  AVTINSP 

441  FORMAT('  AVERAGE  TIME  INSP  =  \F12.4) 

WRITE(*,442)  AVTEXP 

WRITE(4,442)  AVTEXP 

442  FORMATC  AVERAGE  TIME  EXPIR.  =  \F12.4) 

WRITE(*,443)  A  VTTOTAL 

WRITE(4,443)  A  VTTOTAL 

443  FORMATC  AVERAGE  TOTAL  TIME  =  ,F12.4/) 
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WRITER *,444)  VTTIMEINSP 
WRTTE(4,444)  VTTIMEINSP 

444  FORMATC  VT/Ti  -  ',FtL4) 

WRITE(*,446)  TINTOTAL 
WRITE(4,446)  TINTOTAL 
WRITER, 4400)  RFREQ 
WRITE(4,4400)  RFREQ 

4400  FORMATC  BREATHING  FREQUENCY  =  \F12.4) 
WRITE(*,445)  MINVENT 
WRJTE(4,445)  MINVENT 

445  FORMATC  MINUTE  VENT.  =  *,F12.4/> 

446  FORMATC  Ti/Ttot  =  ,tF12.4/) 

WRITE(*,447)  AVINSPPOST 
WRTTE(4,447)  AVINSPPOST 

447  FORMATC  AVERAGE  POST-INSPIRATORY  PAUSE  =  ’.F12.4) 
WRITE(*,448)  AVEXPPOST 

WRITE(4,448)  AVEXPPOST 

448  FORMATC  AVERAGE  POST-EXPIRATORY  PAUSE  =  \F12.4) 
c 

CLOSE(4) 

GOTO  9000 


cccaxcccccccccaxccccccccccccccccccccccccccccccccccccccccccccccccczcccxcccccuxccccx&xcaxcccccc 

c 

c  COMPUTING:  FLOW- VOLUME  RELATIONSHIPS, 
c 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 

6000  OPEN(4,FTLE  =  FTLE2,STATUS  *  ’NEW) 

TOPVOL=0.0 

c 

DO  6  P  =  1,1 
c 

DO  60  J=  1,200 
RVOUN(J)  =  0.0 
VOLEX(J)=0.0 
RINFLOW(J) = 0.0 
RVOLEX(J)  *0.0 
REXFLOW(J) =0.0 
60  CONTINUE 
c 

WRITE(*,63) 
c  WRITE(4,63) 

DO  610  PPP  =  IONSET(P),(EONSET(P)-l) 

PP  *  (PPP-IONSET(P))  + 1 

INVOL = PNEUMO(PPP)-B  ASEFLO  W 

RINVOL  =  FLOAT(INVOL) 

RINVOL = RINVOL/VOLA  VE*(PBAR  +  PREG)/(PBAR)  *310/ 

•  (273+ TA)*(PBAR  +  PLC-PH20)/(PBAR  +  PLC-64.056) 

IF  (RINVOL.LT.(0.0))  GOTO  6100 
RINVOL  *  0.0- RINVOL 
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6100  RVOUN(PP)»RVOLIN(PP-l)-RINVOL 

RINFLOW(PP)  =  FLOAT(lNVOL) 

RINFLOW(PP) = RINFLOW(PP)/VOLAVE*50.0*(PBAR  +  PREG)/(PBAR) 

*  *310/(273 + TA)*(PBAR  +  PLC-PH20)/(PBAR  +  PLC-64.056) 
WRITE(*,631)  PPP,R VOLIN (PP)JUNFLOW(PP) 

WRITE(4,631)  PPP,RVOLIN(PP),RINFLOW(PP) 

610  CONTINUE 
c 

TOP  VOL = RVOLIN(PP) 
c 

WRITE(*,630) 
c  WRITE(4,630) 

c 

DO  620  PPP=STARTEXP(P),(ENDEXP(P)-1) 

PP = (PPP-STARTEXP(P))  + 1 
EXVOL=  PNEUMO(PPP)-BASEFLOW 
REXVOL = FLOAT(EXVOL) 

REXVOL=  REXVOL/VOLAVE*(PBAR  +  PREG)/(PBAR)*310/ 

*  (273 + TA2)*(PBAR  +  PLC-PH202)/(PBAR  +  PLC-64.056) 

RVOLEX(PP)  =  RVOLEX(PP-l)  +  REXVOL 

VOLEX(PP) = TOPVOL-R  VOLEX(PP) 

REXFLOW(PP)  =  FLOAT  (EXVOL) 

REXFLO  W  (PP)  =  REXFLOW(PP)/VOLAVE*50.0*(PBAR  +  PREG)/(PBAR) 

*  *310/(273 +TA2)*(PBAR  +  PLC  PH202)/(PBAR  +  PLC-64.056) 
WRrrE(*,631)  PPP,VOLEX(PP),REXFLOW(PP) 

WRITE(4,631)  PPP,VOLEX(PP),REXFLOW(PP) 

620  CONTINUE 
c 

6  CONTINUE 

c 

63  FORMATO  INSPIRATORY  DATA/ 

*  '  #  Vol(BTPS)  How  ') 

630  FORMAT(/'  EXPIRATORY  DATA  / 

*  ’  #  VoI(BTPS)  How  ) 

631  FORMAT*  150(I4,1X,F8.4,2X,F9.4)) 

CLOSE(4) 

c 

c  Subroutines  for  computing  FEV1.0. 
c 

EVOL-O 
FEV1=0.0 
EXVOLINSP  =  0.0 
C 

DO  64  P  =  STARTEXP(1),(STARTEXP(  1 )  +  49) 

EVOL = PNEUMO(P)-BASEFLOW 
ER  VOL = FLOAT(E  VOL) 

EXVOLINSP  =  EXVOLINSP  +  ER  VOL 

64  CONTINUE 
c 

FEV1  =  EXVOLINSP/VOLAVE*(PBAR  +  PREG)/(PBAR)*310 

*  /(273+TA2)*(PBAR  +  PLC-PH202)/(PBAR  +  PLC-64.056) 
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WRITE(*,65)  FEV1 

65  FORMAT (///■  FEV1  =  ’.PBA,’  Litres  BTPS.’///) 

c 

ceectTecerceTccetrretrrrrrjcrretTexTecrxTX-cxTtTXTececccrcccccccccm-rrr-m-crcrTiTccrrccrrrccc^ccc 

9000  WRITE(*,9996)  FILE1 

9996  FORMAT(/'  ANALYSIS  OF  'tA20,'  IS  COMPLETE.  /) 

WRITE(*,9997) 

9997  FORMATC  WHAT  DO  YOU  WISH  TO  DO  NEXT?’/ 

•  '  1  RETURN  =  ANALYSE  ANOTHER  FILE  WITH  SAME  CALIBRATION,'/ 

*  *  9  RETURN  =  END  Program.’/) 

READ(*,100)  REPEAT 

IF  (REPEATJEQ.l)  GOTO  DIO 
IF  (REPEAT .EQ. 9)  CONTINUE 
c 

9999  STOP 
END 
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APPENDIX  8 


Example  analysis  program  for  use  during  an  exercise  protocol. 


Many  of  the  comments  necessary  to  understand  this  program  and  its  variables  have  been  omitted  where  they 
have  been  provided  in  preceeding  Appendices. 

CCCCCCCCCC<T<TCCwrrfrreeeeffffrffrfrfrrrfrrffrrrffrrfrrcrrrccrrrrrcfrrrrrrrrrrrrrfrerjrrerrefrer 

C 

c  Program  TO  ANALYSE  RESPIRATORY  VOLUME  AND  MOUTH  PRESSURE  DATA, 
c 

c  Program  uses  data  obtained  from  A/D  program 

c  using  the  DASH16  A/D  board  (MetraByte,  at  50Hz), 
c  and  produces  breath-by-breath  and  20  sec  average  output, 

c 

c  ASCII  INPUT  DATA  REQUIRED: 
c  ASCII  file  contains  2  columns  of  data: 

c  Column  1  «  mouth  pressure,  Column  2  =  flow, 

c 

c  THE  FOLLOWING  ANALYSES  ARE  PERFORMED: 
c  (1)  Tidal  volumes,  (2)  Respiratory  timing  (Ti,  Te,  Ttot, 
c  Vt/Ti,  Ti/Ttot,  fb),  (3)  Minute  ventilation,  (4)  Peak-to-peak 
c  mouth  pressures  (cmH20),  and  (5)  flow-volume-pressure  data 
c  for  the  first  breath  of  each  20  sec  interval, 
c 

c  Limitations:  (1)  Expiratory  data  cannot  be  measured  using  this 
c  program  without  changes.  It  has  been  designed 

c  for  a  set-up  using  inspiratory  pneumotach  only, 

c  (2)  Data  is  read  in  from  an  ASCII  file  input  only, 

c  (3)  Convention  for  flow:  inspiration  is  negative, 

c 

c  NIGEL  AJS.  TAYLOR.  June  1990 

c  Diving  Medicine,  Naval  Medical  Research  Institute 

c 

c  LATEST  UPDATE:  10/12/90 
c 

c  This  program  allows  for  unequal  inspired/expired  volumes 
c  by  using  only  complete  cycles,  starting  with  an  inspiration, 

c  Inspiration  is  detected  from  a  negative  signal  from  the  flow 

c  channel.  The  phase  shift  between  the  flow  and  the  pressure 
c  signals  is  taken  into  account  during  analysis, 
c 

c  Three  types  of  data  output  files  are  used: 

c  (1)  file  holding  breath-by-breath  data  (variables  1  to  4); 
c  (2)  file  holding  flow-volume-pressure  data  (20  msec  intervals); 
c  (3)  a  database  holding  only  mean  values  from  variables  in  file  1. 
c  These  correspond  with  Units  7,  6  and  3  respectively  in  program, 
c 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 

$STORAGE:2 
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I 


i 


I 


INTEGER  ZEROFLOWJJTER,BASEUM,PHASEl,PHASE2,INVOL,ZZ2L,EPRESS 
INTEGER  PNEUMO,NOFLOW, FLOW, BAS  E4DATA,COMP,EXVOL,CODE,IPRESS 
INTEGER  GAIN,IONSETEONSET,SEC15,POINT,POINT2,EXPER,LIMBASE 
INTEGER  MAX2,IVOL,BASEFLOW,TBASELIM,EX,EXP,XX,PLOOP,MLOOP 
INTEGER  AACJ),EJsItJ,K4,M,N,P,PP,PPP,Q,R,S,T,W,X,Y^ 

INTEGER  ARBITLIM.PRESS, COUNT, ECOUNT,A2,TIME,LIMITVOL,ELIMBASE 
c 

REAL  VOLUME, V0LAVE,RV0LSUM,PBAR,PH20,TA-AVTIDALV,T1NSP,ETI ME 
REAL  TOTAL VOL,INSPVOL,SYRINGE,VOL,TIDALV,TEXP,TTOTAL/'VRATIO 
REAL  VOLINSP,CALVOURVOL,VTTIMEINSP,AVTINSP,REGCOEFF, CONSTANT 
REAL  AVTEXP,AVTTOTAL,RFREQ,MINVENT,MINV2,TINTOTAL,RSEC15,ZERO 
REALSMALLTV,TIDE,REXVOL,VOLEX3XT5,MAXv\VEXPPOST,TOPVOL,RVOLEX 
REAL  RINVOL,R  VOLIN, RINFLOW,REXFLOW,PINSP,PEXP, INT, EXT, INT5 
REAL  A3,A4,A5tA6,A7,A8^9,A10,All,A12^13, RATIO, EBIT, BIT 
REALM  PEAK,EPEAK,MPPEAK,PPEAK,MPRESS, EMPRESS 
c 

DIMENSION  IDATA(2,1000),LITER(2000),PRESS(1000),MPRESS(900) 

DIMENSION  PNEUMO(1000),VOL(15),VOLBTPS(I5),IONSET(15),TINSP(15) 
DIMENSION  CALVOL(9),EONSET(15))TIDALV(15),TEXP(  15),TTOTAL(15) 
DIMENSION  PINSP(15),PEXP(15),EMPRESS(900) 

DIMENSION  RVOUN(250),RINFLOW(250),RATIO(15) 

DIMENSION  PEAK(900),EPEAK(900)  ,PPEAK(  15) 
c 

CHARACTER  *  10  FILENAME, F1LE1,FTLE2,FILE4,FILE5,FILE6,FTLE7 
CHARACTER  *12  FILE3,ID,A1 

c 

cccccccc<xccccccccccccccccccccccccccccccccccccccaxccccccccccccccccccccccccccccccccccccccccccccccc 

c 

c  INTEGRATION  OF  FLOW  TO  GET  TIDAL  VOLUMES, 
c 

OPEN(3) 

WRITE(*,1) 

1  FORMAT  (///////////////// 

*  GERMAN  1500  fsw  DIVE  SERIES:'// 

*  RESPIRATORY  VOLUME  AND  MOUTH  PRESSURE'/ 

*  '  ANALYSIS  SOFTWARE.'/// 

*  •  win/) 

c 

FLOW  *-200 
SYRINGE  =  000 

ARBITLIM  =  4000 
FILE3  =  GERMBASE.PRN 


100  FORMAT(Il) 

1006  WRITE(MOOO) 

1000  FORMAT(//'  ENTER  NAME  OF  ZERO  FLOW  CALIBRATION  FILE:'//) 
READ(*/20)  FILE1 

WRITE(*,1001) 

1001  FORMAT(//'  ENTER  NAME  OF  HIGH  FLOW  CALIBRATION  FILE:'//) 
READ(*,20)  FILE2 
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WRTTE(*,1002) 

1002  FORMAT(//'  ENTER  NAME  OF  PRESSURE  CALIBRATION  FILE:'//) 
READ(*,20)  FILE4 
C 

C  Reading  calibration  data  from  files. 

C 

OPEN  (6,FTLE = FILE4, STATUS  =  'OLD') 

READ(6,50)  REGCOEFF, CONSTANT 
CLOSE(6) 

WRTTE(*,5099)  REGCOEFF, CONSTANT 
5099  FORMAT//1  REGRESSION  COEFFICIENT  (cmH20)  =  ’,F92/ 

•  1  CONSTANT  A  (cmH20)  =  ',F9 2/f) 

c 

OPEN(3,FILE= FILE1, STATUS  =  OLD) 

READ(3,101)  ZEROFLOW 
CLOSE(3) 
c 

WRITE/ *,102)  ZEROFLOW 

101  FORMAT(I7) 

102  FORMATC  BASAL  FLOW  (RAW  UNITS)  «  ’,15//) 

50  FOR  MAT(2X,F9.2^2X,F9.2) 

c 

BASE  =  (FLOW  -  ZEROFLOW)  *4/100 
BASELIM  =  BASE  +  ZEROFLOW 
TBASEUM  *=  ZEROFLOW  -  BASE 
C  WRITE(*,103)  BASE 

103  FORMATC  VALUE  OF  BASE  =  15//) 

C  WRITE/*, 104)  BASELIM 

104  FORMATC  VALUE  OF  BASELIMIT  *  15//) 

C 

OPEN(3,FTLE = FILE2.STATUS  =  OLD',FORM  =  UNFORMATTED) 

READ(3)  (UTER(jy  =  L2000) 

11  FORMAT(10(I7)) 

CLOSE(3) 

c 

M  =  1 
L=0 

12  DO  120  J  =  M.2000 

IF  (UTER(J).LT .BASELIM)  GOTO  1200 

120  CONTINUE 

1200  CONTINUE 

IF  (J.GE2000)  GOTO  129 

c  Integration  of  flow  data. 

VOLINSP=0.0 

121  DO  127  T-J.2000 

IVOL= LITER  (T)-ZEROFLOW 
RVOL=FLOAT(TVOL) 

VOLJNSP  -  VOUNSP-R  VOL 
IF  (LITER (T).GE.BASELIM)  GOTO  1270 
127  CONTINUE 

1270  CONTINUE 

c  Checking  data  for  false  volumes  before  storing  volumes. 
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IF  (VOLINSP-LE-ARBITLIM)  GOTO  128 
L=L+1 

CALVOL(L)  =  VOLINSP 
C  WRITER, 1391)  L,CALVOL(L) 

1391  FORMATC  FOR  PUMP  43,’  CALVOL  =  ’.F12.1) 

128  IF  (T.GE.2000)  GOTO  129 
M  =  T 

GOTO  12 

c  Calculation  of  average  calibration  value. 

129  RVOLSUM=0.0 
DO  1280  W=1,L 

RVOLSUM  =  RVOLSUM  +  CALVOL(W) 

1280  CONTINUE 

VOLAVE  =  RVOLSUM  /  (SYRINGE *L) 

1281  FORMAT(10(I7)) 

WR1TE(»,1282)  VOLAVE 

1282  FORMAT(/’  MEAN  FLOW  CALIBRATION  VALUE  =’,F14.4,'  bits/litre’//) 
WRITE(*,1360) 

1360  FORMAT(//’  LOOK  AT  THE  VOLUME  CALIBRATION  VALUE:’/ 

*  ’  SELECT  A  FALSE  VOLUME  LIMIT  BASED  UPON  THIS  VALUE:’/ 

*  *  [e.g.  200  ml  =  20%]’// 

*  ’  ENTER  YOUR  VALUE  {17  -  NO  decimal}’//) 

READ(*,1361)  LIMITVOL 

1361  FORMAT(I9) 

C 

2  FORMAT(//’  ENTER  NAME  OF  DATA  FILE  TO  BE  ANALYSED:’///) 

20  FORMAT(AIO) 

14  FORMAT(///’  HOW  MANY  MINUTES  DID  THIS  EXPERIMENT  RUN?  [13]’//) 
WRITE(*,130) 

130  FORMAT(/’  WHAT  WAS  THE  BAROMETRIC  PRESSURE  (mmHg  F12.2)?’/) 
READ(*,131)  PBAR 

131  FORMAT  (F122) 

PBAR  =  PBAR  *  136 
WRITER, 132) 

132  FORMAT(//’  WHAT  WAS  THE  INSPIRED  AIR  TEMPERATURE  (degC  F8.4)?’/) 
READ(*,131)  TA 

WRITE(*,133) 

133  FORMAT(//’  WHAT  WAS  INSPIRED  WATER  VAPOUR  PRESSURE  (mmHg  F8.4)?’/) 
READ(*,131)  PH20 

PH20  =  PH20  *  136 


COMP  =  0 
COMPRESS  =  0.0 
PLC  =  0 
SIPLC  =  0 
PREG  =  0 
PSTERN  -  0 
PHASE1  =  0 

c  PHASE1  comment: 

c  The  flow  and  mouth  pressure  signals  are  rarely  in  phase  so 
c  this  difference  must  be  taken  into  account.  Let  the  value  for 


the  flow  signal  be  zero,  then  subtract  the  number  of  points 
between  the  2  signals  from  zero.  Report  the  sign  below. 


c 
c 
c 

WRITER, 8905) 

8905  FORMAT(//'  WHAT  WAS  THE  PHASE  SHIFT  BETWEEN  FLOW  AND  MOUTH' 

*  '  PRESSURE  [FLOW  MINUS  PRESSURE:  I3]?'//) 

READ(*,135)  PHASE1 
C 

135  FORMAT(B) 

1444  FORMAT(I2) 

1310  WRITE(*,2) 

READ(*,20)  FILENAME 
WRITE(M4) 

READ(*,1444)  SEC15 
WRTTE(,,8901) 

8901  FORMAT(//'  ENTER  SUBJECTS  ID  CODE:'//) 

READ(*,20)  ID 

WRITE(*,8904) 

8904  FORMAT(//'  ENTER  CONDITION  CODE  (numerical  13):'//) 

READ(*,101)  CODE 
WRITE(*,8902) 

8902  FORMAT(//'  ENTER  NAME  OF  FILE  FOR  FLOW-VOL-PRESSURE  OUTPUT:'//) 
READ(*,20)  FILE6 

WRITE(  *,8903) 

8903  FORMAT(//'  ENTER  NAME  OF  FILE  FOR  BREATH-BY-BREATH  OUTPUT:'//) 
READ(*,20)  FILE7 

C 

c  Reading  formatted  file. 

c  Placing  flow  data  into  the  array  PNEUMO  and  pressure  data  into  PRESS 

C 

OPEN(3,FILE  =  FILENAMESTATUS  =  OLD, FORM  =  FORMATTED) 

OPEN(5,FILE  =  FILE3, STATUS  =  OLD'.FORM  =  'FORMATTED') 

DO  1592  B  =  1,10000 

READ(5,L590,END=  1591)  Al,A2A3,A4,A5,A6,A7,A8,A9rA10,All,A12y4.13 
1592  CONTINUE 
1591  CONTINUE 
1590  FORMAT  (A6,I5, 1 1F83) 

ENDFILE  5 
BACKSPACE  5 

OPEN(6,FILE  -  FILE6, STATUS  =  NEW, FORM  =  FORMATTED) 

OPEN(7,FILE  =  FILE7, STATUS  =  'NEW, FORM  =  FORMATTED) 
c 

TIME  =  0 
ETIME  *  0.0 
PLOOP  =  0 
SEC15  =  SEC15  *  3 
MLOOP  *  SEC15 
c 

CCCCCCCCCCCCCCCCCCfCCCCCfCCCCCCCCCCCCCCCCCCCCCCCCCCCCJCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 

DO  4446  PLOOP  -  1, MLOOP 
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r 


c 


c 


21 

2090 

2999 

2998 

C 


22 

C 

C 

2292 
C 

2293 
C 

C 

24 

• 

c 

c 

6003 


ETIME  =  ETIME  +  03334 

TIME  =  TIME  +  1 

POINT  =  1 

POINT2  =  1000 

BASEFLOW  =  ZEROFLOW 

J=0 

K=0 


DO  21  K = POINT, POINT2 
READ(3,2090,END  =  2999)  (IDATA(L,K),L=  1,2) 

CONTINUE 
GOTO  2998 
FORMAT(2I6) 

WRITEC34)  FILENAME 
CONTINUE 

DO  22  J  =  POINT, POINT2 
PRESS(J)  =  IDATA(U) 

PNEUMO(J)  =  IDATA(2rI) 

CONTINUE 

WRITEC3292)  IDATA(1, POINT) 

FORMATO  STARTING  POINT  FOR  20  SEC  ,15) 
WRITE(*,2293)  J 

FORMATC  FINISHING  POINT  FOR  20  SEC  ,14) 

GOTO  6003 

FORMATC  END  OF  FILE  ENCOUNTERED.1/ 

ANALYSIS  OF  RLE:  >10,'  COMPLETE.1//) 


Volume  analysis  routine. 

N  =  0 


c 

DO  3002  J  =  1,15 
TIDALV(J)  =  0.0 
VOLBTPS(J) = 0.0 
IONSET(J)  =  0 
EONSET(J)=0 
TINSP(J)=0.0 
TEXP(J)=0.0 
TTOTAL(J)  =  0.0 
3002  CONTINUE 

TOTAL  VOL  =  0.0 
c 

c  COUNTERS:  I  =  inspiration  onset,  E  =  expiration  onset, 
c  P  =  pairs  of  insp  and  exp,  J,K,M  =  starting  points, 

c  L  *  tidal  volumes. 

c  NOTE:  I  is  used  as  a  respiratory  cycle  counter  from  this  point. 
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c  Integrating  flow  to  obtain  volume  (inspiration  -ive). 
c 

E=0 

1=0 

J=0 

K=0 

L=0 

M  =  1 

C 

3  CONTINUE 

C 

DO  3111  J  =  M.POINT2 
IF  (PNEUMO(J).GE.BASELIM)  GOTO  3112 
C  WRITE(*,8820)  PNEUMO(J),J 

3111  CONTINUE 

88220  FORMAT('+  PNEUMO(J)  =  ',16,'  FOR  J  =  ,16) 
C  WRITE(*,8822) 

8822  FORMAT(//'  FINISHED  LOOP '//) 

C 

3112  DO  31  K=J,POINT2 

IF  (PNEUMO(K).LT.BASELIM)  GOTO  310 
C  WRITE(*,8820)  PNEUMO(J),J 

31  CONTINUE 

C 

310  CONTINUE 

C  WRITE(*,8822) 

IF  (K.GE.POINT2)  GOTO  3102 
VOLINSP  =  0.0 
EXVOLINSP =0.0 


DO  32  T  =  K,POINT2 

IF  (PNEUMO(T)  .GE.B ASELI M )  GOTO  33 
I  VOL = PNEUMO(T)-BASEFLO  W 
RVOL=FLOAT(IVOL) 

VOLINSP  =  VOLINSP  -  RVOL 

32  CONTINUE 

33  CONTINUE 

c  The  -RVOL  converts  volume  *.o  positive, 

c  Check  for  small,  false  volumes, 

c 

IF  (VOLINSP.LT ARBITLIM)  GOTO  3103 
c 

c  Marking  points  of  inspiration  start  (IONSET),  &  end  (EONSET), 
c  and  expiration  start  (STARTEXP),  &  end  (ENDEXP). 
c 

1  =  1  +  1 

IONSET(I)  =  K 
EONSET(I)  =  T 
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c 


330 

DO  330  J  =  1,15 

VOL(J)=0.0 

CONTINUE 

c 

L=L+1 

3390 

FORMATC 

VOLUME  INSPIRED 

VOL(L)  =  VOLINSP  /  VOLAVE 


C  WRITE(*3390)  VOL(L) 
c 

c  Volume  BTPS  correction  (deg  C  and  cmH20  units) 

c 

VOL(L)  =  VOL(L)  *  (PBAR  +  PREG)  /  (PBAR) 

VOLBTPS(L)  =  VOL(L)  *310 /(273  +  TA)*(PBAR  +  PLC-PH20) 

*  /  (PBAR  +  PLC-64.056) 

TIDAL V(L)  =  VOLBTPS(L) 

4199  FORMAT  (316, F83) 

c 

3103  IF  (T.GE.POINT2)  GOTO  3102 

C  WRITE(*,4199)  L,IONSET(L),EONSET(L), TIDAL  V(L) 

IF  (VOLINSP.GTARBITLIM)  GOTO  3104 
M  =  T 
GOTO  3 
c 

3104  CONTINUE 
E  =  E+1 

3410  CONTINUE 
M=T 
GOTO  3 

3102  CONTINUE 
c 

8810  FORMATC  HERE',15) 

8000  FORMATC  ANALYSIS  OUTPUT  FROM  DATA  FILE:  ,A15/) 
4  FORMATC  START/INSP  END/INSP  TIDAL  VOLUME  ) 

41  FORMAT(3I63X,2F12.4) 


c 

c  RESPIRATORY  TIMING  SUBROUTINES  (including  post-insp 

c  [POSTTNSP]  and  post-expiratory  pauses  [POSTEXP]). 

c 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 


c 


c 


IF  (I.GT.l)  ZZZ  =  (1-1) 
IF  (I.LE.l)  ZZZ  =  1 


DO  42  P=1,ZZZ 

TINSP(P)  =  (EONSET(P)  -  IONSET(P))  /  50.0 
IF((I-1).LE.0)  GOTO  4271 
TTOTAL(P)  =  (IONSET(P+  1)-IONSET(P))/50.0 
IF  (TTOTAL(P).LE.(O.O))  GOTO  4271 
GOTO  4272 
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4271  TTOTAL(P)  =  TINSP(P)  *  2 

4272  TEXP(P)  =  TTOTAL(P)  -  TINSP(P) 

42  CONTINUE 

C 

422  FORMAT(15(F7.4,5X,F7.4,lX,F7.4)) 

C 

AVTINSP =0.0 
AVTEXP=0.0 
AVTTOTAL =0.0 
AVTIDALV  =  0.0 
c 

DO  44  P  =  1,ZZZ 

IF  ((I-l).LE.O)  GOTO  4007 

AVTIDALV  =  AVTIDALV  +  TIDALV(P) 

AVTTNSP  =  AVTINSP  +  TINSP(P) 

AVTEXP  =  AVTEXP  +  TEXP(P) 

AVTTOTAL  =  AVTTOTAL  +  TTOTAL(P) 

44  CONTINUE 
C 

AVTIDALV  =  AVTIDALV  /  (1-1) 

AVTINSP  =  AVTINSP  /  (1-1) 

AVTEXP  =  AVTEXP  /  (1-1) 

AVTTOTAL  =  AVTTOTAL  /  (1-1) 

IF  (AVTTOTAL.LE.(O.O).OR-AVTINSP.LE.(O.O))  GOTO  4009 
RFREQ = 60.0/ AVTTOTAL 
VTTIMElNSP  =  AVTIDALV  /  AVTINSP 
MINVENT  =  AVTIDALV  *  RFREQ 
TINTOTAL  =  AVTINSP  /  AVTTOTAL 
GOTO  4008 

4007  AVTIDALV  =  TIDALV(l) 

AVTINSP  =  TINSP(l) 

AVTEXP  =  TEXP(l) 

AVTTOTAL  =  TTOTAL(l) 

IF  (TTOTAL(1).LE.(0.0).OR.TINSP(1).LE.(0.0))  GOTO  4009 
RFREQ  =  60.0  /  TTOTAL(l) 

VTTIMElNSP  =  TIDALV(l)  /  TINSP(l) 

MINVENT  =  TIDALV(l)  *  RFREQ 
TINTOTAL  =  TINSP(l)  /  TTOTAL(l) 

GOTO  4008 

4009  RFREQ  *  RFREQ 

VTITMEINP  =  VTTIMElNSP 
MINVENT  =  TIDALV(l)  *  P  -  TEQ 
TINTOTAL  =  TINTOTAL 

4008  CONTINUE 
c 

ccccccccccccccccccc<xcccccccccccccaxcccccccccx'cccccx~x'xrccccccccccccc«xxcccccccccccccccccccccccccc 

c 

c  FLOW- VOLUME-PRESSURE  SUB-ROUTINES 

c  Inspiratory  output  only,  for  the  first  and  last  breaths  and  for 
c  the  First  breath  in  the  middle  20  sec  period  of  each  2  minutes. 
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TOPVOL=0.0 
BIT  =  0.0 
EBIT  =  0.0 


DO  500  P  =  1,ZZZ 

IF  (P.NE.l)  GOTO  500 

DO  501  J  =  1,250 
RVOLIN(J)=0.0 
RINFLOW(J)  =  0.0 
CONTINUE 

DO  5011  J  =  1,900 
MPRESS(J)=0.0 
PPEAK(J)=0.0 
EMPRESS(J)  =0.0 
EPEAK(J)  =  0.0 
CONTINUE 


DO  510  PPP=10NSET(P),(EONSET(P)-l) 

PP = (PPP-IONSET(P))  + 1 

BIT  =  BIT  +  0.02 

IPRESS  =  PRESS(PPP + PHASE1) 

MPRESS(PP)  =  FLOAT(IPRESS) 

MPRESS(PP)  =  (MPRESS(PP)-CONSTANT)/REGCOEFF 
INVOL=  PNEUMO(PPP)-BASEFLOW 
RINVOL=  FLOAT(INVOL) 

RINVOL=  RINVOL/VOLAVE*(PBAR  +  PREG)/(PBAR)*310/ 

(273  +  TA)*(PBAR  +  PLC-PH20)/(PBAR  +  PLC-64.056) 

IF  (RINVOL.LT.(O.O))  GOTO  5100 
RINVOL = O.O-RINVOL 
RVOUN(PP)  =  RVOLIN(PP-l)-RINVOL 
RINFLOW(PP)  =  FLOAT(INVOL) 

RINFLOW(PP)  =  RINFLOW(PP)/VOLAVE*50.0*(PBAR  +  PREG)/(PBAR) 
•310/(273  +  TA)*(PBAR  +  PLC-PH20)/(PBAR  +  PLC-64.056) 

IF  (TIME.EQ.1)  GOTO  5110 
IF  (TIME.EQJ5)  GOTO  5110 
IF  (TIME.EQ.il)  GOTO  5110 
IF  (TIME.EQ.17)  GOTO  5110 
IF  (TIME.EQ23)  GOTO  5110 
IF  (TIME.EQ29)  GOTO  5110 
IF  (TIME.EQ35)  GOTO  5110 
IF  (TIME.EQ.41)  GOTO  5110 
IF  (TIME.EQ.47)  GOTO  5110 
IF  (TIME.EQ^3)  GOTO  5110 
IF  (TIME.EQJ9)  GOTO  5110 


IF  (TIME.EQ.65)  GOTO  5110 
IF  (TIME.EQ.71)  GOTO  5110 
IF  (TTME_EQ.75)  GOTO  5110 
GOTO  5111 

5110  CONTINUE 

c  WRITE(*431)  BIT,RVOLIN(PP),MPRESS(PP),RINFLOW(PP) 

WRITE(6,531)  BIT,RVOLIN(PP),MPRESS(PP),RINFLOW(PP) 

5111  CONTINUE 

510  CONTINUE 

c 

c  TOPVOL=  RVOLIN(PP) 

c 

c  DO  520  PPP = STARTEXP(P),(ENDEXP(P)-1) 

c  PP  =  (PPP-STARTEXP(P))  + 1 

c  EBIT  =  EBIT  +  0.02 

c  IPRESS  =  PRESS(PPP + PHASE1) 

c  EMPRESS(PP)  =  FLOAT  (IPRESS) 

c  EMPRESS(PP)  =  (EMPRESS(PP)-CONSTANT)/REGCOEFF 

c  EXVOL=PNEUMO(PPP)-BASEFLOW 

c  REXVOL = FLOAT(EXV  OL) 

c  REXVOL  =  REXVOL/VOLAVE*(PBAR  +  PREG)/(PBAR)*310/ 

c  *  (273 + TA2)*(PBAR  +  PLC-PH202) /(PBAR  +  PLC-64.056) 

c  R  VOLEX(PP)  =  RVOLEX(PP-l)  +  REXVOL 

c  VOLEX(PP) = TOP  VOL-R  VOLEX(PP) 

c  REXFLOW(PP)  =  REXFLOW(PP) /VOLAVE*50.0*(PBAR  +  PREG)/(PBAR) 

c  *  *310/(273  +  TA2)*(PBAR  +  PLC-PH202)/(PBAR  +  PLC-64.056) 

c  WRITE(*,531)  EBIT1VOLEX(PP),EMPRESS(PP),REXFLOW(PP) 

c  WRITE(6431)  EBrr,VOLEX(PP),EMPRESS(PP),REXFLOW(PP) 

520  CONTINUE 

c 

500  CONTINUE 

531  FORMAT  (F83r3X,F83^X,F83,2X,F93) 


MOUTH  PRESSURE  CALCULATIONS 

CCCCCCCCOXCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCl 

COUNT  =  1 
ECOUNT  =  1 
MPPEAK=0.0 
AVRATIO =0.0 
ZERO  =  CONSTANT 

DO  60  J  =  l,900 

MPRESS(J)=0.0 
EMPRESS(J)  =0.0 
PEAK(J)  =  0.0 
EPEAK(J)=0.0 
60  CONTINUE 
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mu 


DO  6111  J  =  l,15 

PPEAK(J)=0.0 
RATIO(J)=0.0 
6111  CONTINUE 

c 

DO  6  P  =  1,ZZZ 
c 

DO  600  PPP = IONSET (P),(EONSET (P)-l) 

PP  =  (PPP-IONSET(P))  + 1 
IPRESS = PRESS  (PPP  +  PHASEl) 

MPRESS(PP)  =  FLOAT  (IPRESS) 

MPRESS(PP)  =  (MPRESS(PP)-CONSTANT) /REGCOEFF 
IF  (MPRESS(PP).GE.PEAK(COUNT))  GOTO  6000 
COUNT  =  COUNT  +  1 
PEAK(COUNT)  =  MPRESS(PP) 
c  WRITER, 6401)  PEAK(COUNT) 

6000  CONTINUE 

600  CONTINUE 

c 

c  IF  (PEAK(COUNT).LT.(-IOO.O))  PEAK(COUNT) = 0.0 

c  WRITER, 6400)  PEAK(COUNT) 

FORMAT('  PEAK  INSPIRATORY  PRESSURE  =  ,F93) 

FORMATC  PEAK  EXPIRATORY  PRESSURE  =  \F93) 

FORMATC  +  INSPIRATORY  PRESSURE  REAL  =  ',F9.1) 
FORMATC  +  EXPIRATORY  PRESSURE  REAL  =  \P9.1) 
FORMAT('  PEAK-TO-PEAK  =  F9.1'  MEAN  =  ’F9.1'  BREATH  ,13/) 

DO  610  PPP=EONSET(P),(IONSET(P+ 1)-1) 

PP  =  (PPP-EONSET  (P))  + 1 
EPRESS  =  PRESS(PPP  +  PHASEl) 

EMPRESS  (PP)  =  FLOAT  (EPRESS) 

EMPRESS(PP)  =  (EMPRESS(PP)-CONSTANT) /REGCOEFF 
IF  (EMPRESS(PP).LE.EPEAK(ECOUNT))  GOTO  6100 
ECOUNT  =  ECOUNT  +  1 
EPEAK(ECOUNT)  =  EMPRESS(PP) 

WRITE(*,6403)  EPEAK(ECOUNT) 

CONTINUE 
CONTINUE 

IF  (EPEAK(ECOUNT).GT.(IOO.O))  EPEAK(ECOUNT)=0.0 
WRITE(*,6404)  EPEAK(ECOUNT) 

PPEAK(P)  =  EPEAK(ECOUNT)-PEAK(COUNT) 

WRrTE(*,6402)  PPEAK(P),P 

FORMATC  PEAK-TO-PEAK  PRESSURE  =  \F93,'  BREATH  I4) 
IF  (EPEAK(ECOUNT).LE.(O.O))  GOTO  6138 
RATIO(P)  =  (O.O-PEAK(COUNT))  /EPEAK(ECOUNT) 

GOTO  6139 
RATIO(P)  =  0.0 
CONTINUE 

WRITE(7,41)  P,IONSET(P),EONSET(P),TIDALV(P),PPEAK(P) 


c 

6100 

610 

c 

c 

c 

c 

c 

6402 


6138 

6139 
c 


68 


c 


6288 


6289 

c 

6 

c 


IF  ((I-l)JLE.O)  GOTO  6288 
MPPEAK  =  MPPEAK  +  PPEAK(P) 
WRITE(*,6405)  PPEAK(P),MPPEAK,P 
AVRATIO  =  AVRATIO  +  RATIO(P) 
GOTO  6289 
CONTINUE 
MPPEAK  =  PPEAK(l) 

AVRATIO  =  RATIO(l) 

CONTINUE 


CONTINUE 


MPPEAK  =  MPPEAK  /  ZZZ 
AVRATIO  =  AVRATIO  /  ZZZ 


IF  (PLOOP.GT.l)  GOTO  4445 

WRITE(  *  ,4444) 

FORMAT(// 

AVERAGE  VALUES  IN  20  SECOND  SEGMENTS’/ 

Vt  P-Peak  Ti  Te  Ttot  Vt/Ti  ’ 

'  lb  Vc  Ti/Ttol’) 

WRITF (*,1599)  A VTIDALV, MPPEAK, 
AVTINSPAVTEXP,AVTTOTAL,VTTIMEINSP,RFREQ, 

MINVENT.TINTOTAL 
FORMAT(9F8J) 

WRITE (5, 1590)  ID,CODE,ETIME,AVTlDALV,MPPEAICAVRATIO, 

•  AVTINSPAVTEXP^VTTOTAL,VTTIMEINSP,RFREO, 

*  MINVENT.TINTOTAL 

c 

c 

4446  CONTINUE 

c 

c 

ccccccccccccccccccccccc  ccmxcccaxccccccccccccccccaxcccccccccccccccccaxcccccccccccccccccccccccc 


4444 


c 

4445 


1599 


9995  WRITE(\999t)  FILENAME 

9996  FORMAT(///  ANALYSIS  OF  FILE  ’.A10,’  IS  COMPLETE.'///) 

CLOSE(3) 

CLOSE(5) 

CLOSE(6) 

CLOSE(7) 

A=0 

WRITE(*,9997) 

9997  FORMAT(’  WHAT  DO  YOU  WISH  TO  DO  NEXT?’/ 

•  ’  l  RETURN  *=  ANALYSE  MORE  DATA  USING  SAME  CAL.  FILES,’/ 

•  ’  6  RETURN  =  ANALYSE  MORE  DATA  USING  DIFFERENT  CAL.  FILES,'/ 

•  ’  9  RETURN  »  TERMINATE  Program.’/) 

READ(*,9998)  A 

9998  FORMAT(Il) 


9999 


IF  (A.EQ.9)  GOTO  9999 
IF  (A£Q.l)  GOTO  DIO 
IF  (A.EQ.6)  GOTO  1006 
STOP 
END 
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PERMS  AND  ABBREVIATIONS 


Lung  volume  and  ventilation  variables  (liters  BTPSV 
ERV  =  expiratory  reserve  volume 

FRC  =  functional  residual  capacity 

1C  =  inspiratory  capacity 

DRV  =  inspiratory  reserve  volume 

RV  =  residual  volume 

TLC  =  total  lung  capacity 

VC  =  vital  capacity 

VT  =  tidal  volume 

ve  =  minute  ventilation  derived  from  expiratory  data 

MVV  =  maximal  voluntary  ventilation 

EEVi  0  =  forced  expired  volume  during  first  second  of  maximal  expiration 

FEVt.o/VC  =  ratio  of  FEVj  0:vital  capacity  expressed  as  a  percentage. 

Respiratory  flow  variables  fliters.second 

=  peak  expiratory  flow  rate 
=  maximum  expiratory  flowrate  at  25%  of  VC 
=  maximum  expiratory  flowrate  at  50%  of  VC 
=  maximum  expiratory  flowrate  at  75%  of  VC. 


PEFR 

MEFRjs 

MEFRm 

MEFRtj 


Respiratory  timing  variables  (seconds!: 

T]  =  time  of  inspiration 

Tg  *  time  of  expiration 

tTOT  -  total  respiratory  duration 
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Post-exp  pause  =  duration  of  the  pause  between  the  end  of  expiration  and  the  commencement  of  the 


Post-insp  pause 

ti/ttot 

vT/r, 


next  inspiration 

=  duration  of  the  pause  between  the  end  of  inspiration  and  the  commencement  of  the 
next  expiration 

=  ratio  of  the  time  of  inspiration  to  the  total  cycle  duration 
=  ratio  of  the  tidal  volume  to  the  inspiratory  duration  (1-s'1). 
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